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Abstract 


An introductory survey of the global energy-balance clisate 
laodels is presented with an cssphasia on analytical results ♦ A sequence 
of increasingly cossplicsted esodela involving ice cap and radiative 
feedback processes are solved and the solutions end paraaster senel- 
tivitiea are etudiad. The esodal parasactoriaationa are eHanined 
critically in light of ®any current uncertainties. A slssple seasonal 
®jd®l is used to atudy the effects of changes in orbital eleaants 
on the temperature field. A linear etafeUlty theorca and a cosspXete 
nonlinear otabllity analysis for the taodels are developed. Analytical 
solutions are also obtained for the linearized nsodels driven by 
atochastic forcing eleoents. In this context the relation between 
natural fluctuation statistics and clitaiite sensitivity is stressed. 



Int*:ed«ction 


Tlte theory of cllmt® ha« received isaich atteatioa ia the lf^*t 
years. %« evolution of high-speed cesputere and the devalopRsnt of 
aussarlcal wether prediction ®odela have mde the 8l«ui®tl«n of 
cllaaee at least ponderable. The global ecale collection aud analysis 
of observations have provided a base for developing aiwS verifying 
Bodela. In addition, the extraction of paleoclitaatic lafonBatioa 
froa vaiious sources is beginning to produce a legible* record of 
cliaatic history. Preliminary studies of the earth’s climate 
suggest that the present state my be a delicate one pcasibly vul- 
nerable to unintentional adverse changes by man’s activities. As it 
demands more and more froa the earth’s dvlndllng resources, the 
grm.ing human population becomes less able to cope with climatic 
change. As a result, the forecasting of future climates becomes 
increasingly important. Evidently, the time has come to develop 


mathematical models of the climate, 

Tlie purpose of this paper is to present an Introductory survey 

of simple climate models based upon elciuentary heat balance con- 
siderations. The paper Is intended to be pedagogical, introducing 
complicated subjects by way of solvable examj^les. The paper begins 
with the fundamental principles that govern planetary climates and 
proceeds to develop a theory for onc-dlmenslonal climate models. 

Schneider and Dickinson <19741 have surveyed many approaches to 
climate «>daling. Foremost in their diacusslon is the hierarchy of 
climate models. A wide range of ««odels may be constructed based upon 





the choice end nc®bar of degress of freedom to be Included. Ceaeral 
circulation i^deXe (€CMa) can have up to a aillion degrees o£ freedoa 
vhiie the Bimplest models have only a few. Although in principle 
GCiU could becosas physically realistic* they are ©sspensive ®«d 
cmberso&ws. Ihe artificial climtes generated by these etodels are 
typically as complicated and inscrutable as the earth’s climte. 

Their major advantage (once perfected) will be their controllability; 
i.e.* the possibility of testing hypotheses by changing boundary 
conditions* a lustury not afforded by the real climate. 

Because of the expense and complicated output of large GCMs* a 
variety of simplified initial value models are presently under 
construction for use in sensitivity expcrlosants (e.g*» Held and 
Suarez* 1978; Gates and Schleslnger* 1977). It is hoped that these 
models will produce climates and climate change responses similar 
to those of their larger counterparts. 

Since the climate is represented by the long-term averages of 
atmospheric variables, one promising approach Is to construct 
equations in tenas of these averaged quantities. Recently* there 
have been many attempts to construct and study such equations. The 
resulting models are referred to as "statistical dynamical taodels" 
(SDMs) and were reviewed by Saltzman (1978). 

Ataong the SDMs are the few variable models. If only the ver- 
tical dimension is retained, one obtains the radiative-convective 
models reviewed by Samnathan and Coaklsy (1978). The advantage of 
these models is that they can be used to compute radiative transfer 



1„ d«»ll and, tharafor., cltotlc teadback B>cb.nt.,a. .».ocl«.d 
with radiation my b® cerefuliy etudled. 

If ». chaiacMicla* • a»l“» «' “ 

.IngU »»=b«T. a«y tka ••« t. ». davalcp ~d.la »ith 

only horUontal dl-nnioM. Zon.1 .v.aglnB anch nad.Ia la.da ua W 
t„. ona-dl-nalonal aU»t. »dal». Modal, of .hi. OT. hav. ba». 
atudlod foi aoffie yaara (Snsatroo, 1928; Fritz. I960; Oplk, 196-', 
EtlkOBOn. 1969). 8ana»ad Intoraat In than, taodal. .a. atlBnlatad by 
Bodyko (1968. 1969. 1972) and Sailor. (1969). 12.... too Invaa.l- 

sator. lndap«.dantiy darWad ona-dlnanslonal K.d.1. ba.ad npon the 
thormodynanlc ..nation. Back van. la th. ..nation vrltt.n In 
tarn, of tb. aan-laval tawaratnr. Ilald. In doing an .hay dlatlllad 
the cllnata problen Into « oaa-dliMnsion«l. steady-stata. bonndarj 
vain, problan that vaa aolvad for the tanparatnta field. M.hongh 
th. aqnatlona vara nonllnaar. aoln.lona aonld be aztr.ctad by al.har 
analytical or rall.bla nnn.rla.1 procadnras. On. «.s than In a 
poaltlon to vary anch -glvan" patanat.r. oa tha aolar aon.tant to 

Study the ssodel response. 

Bodyho and Sailor, arrived at tha aurfaca tanpar.tnra dapandanaa 
of tha Individual tanaa In tha anargy-bolanca mnatlona th.ough 
indapandan. studio, of tha obsatvad heat flnnaa. hot aorprlalngly. 
th. functional fonaa nar. rather different from other, havar- 

thalaaa. both nodal, yielded th. present cllnata as solution.. 
Furthamora. both predicted th. a.na high aao.ltlvlty to change. In 
tha aolar constant. They predicted that If tha .ol.r conava,,- »ar. 



lowered by only & few percent, the polar ice capa would expand 
catastrophically utttil the globe bec®»® coispletely covered by Ice. 

Iba Budyko-Sellers »odels fora the basts for this review, 

A nusHber of questions isaasdiately arise froa the pioneering 
studies just s^sntloned. For instance, to what extent are the aodels 
equivalent to each other and, for that asatter, do they behave like 
sore sophisticated eodels possibly iaitatlng the earth’s climate? To 
what extent is the extreme clioate sensitivity of these models 
dependent upon the paraoeterizatlons used to relate the surface 
temperature to the heat fluxes? What is the nature oi the model 
climate solutions; for instance, is the model solution unique; is it 
stable to small perturbations? Arc the models consistent with the 
history of the solar system? What is the range of space, and time 
scales for which the models are valid? Can the models ^luggest new 
measurements or data reduction methods that would further the develop- 
ment of a climate theory? Are there any ways to test tKse validity of 
the models; for Instance, would they apply to the other planets? Can 
the models be extended to Include seasonal and regional" effects? Can 
the models test various theories of the ice-ages? Finally, can the 
models be used to define research problems for the more -comprehensive 
models? Some of these questions have been answered in tSie last few 
years while others remain open. 

One appealing feature of the Budyko-Sellers models ils their 
simplicity* Ihis simplicity facilitates the use of the jjziodels as 
teaching tools* It# therefore^ seeias appropriate in thisB review 






to use eleMiXitery ©nalytic^l is^tbods t^honawr possiibie* In this vay 
wo can keep tho physical »chanis» feefot'o m at all times* As a 
result » the inevitable fudge factors will be esi^licit* More emphasis 
will b© placed upon theoretical rather than nusserical results becetise 
the latter are subject to change with nm observations and new 
developMuts in para^terigation theory* Althou^ the theory of 
energy-balance models has had mny contributors with varied approaches 
we have atte^apted to develop the theory whenever possible in a simple 
mlforta manner » often drawing upon the work of others. 

The paper is divided into nine flections. After reading Section 2 
most of the other sections can be road independently of each other. 


2* Introduction to heat balatice leodels 

To begin, wa Introduce the concept of global radiative heat 
balance. For slEplicity we asots® that the earth emita radiation 
like e black body, la radiative equilibriusj the rate at which solar 
radiation io ebejorbed asatches the rate at which infrared radiation is 
eisitted. The condition of radiative equllibrlvoa is given by 

47tR^C3T_^ - a (1 “ o„) . (1) 

£1 op 

where is the effective radiating teB 5 >eraturc of the planet; R is 

the radius of the planet; a is the solar constant, taken in this 

o 

-2 

section to be 1340 W a ; o is the Stefaa-Boltsmann constant, 

0.56687 s 10“^ W nT^ K~^; is the planetary albedo defined as 

o - ~ [ dx S(x) a(x) , (2) 

P 2 

where x “ sine of latitude; a(x) Is the albedo r latitude x; and S(x) 
X8 the laean annual distribution of radiation reaching the top of the 
atmosphere normalized so that the integral of S(x) from 0 to 1 is 
unity. S(x) and its seasonal analog S(x,t) may be computed exactly 
(Sellers, 1965), but for this discussion we may use the approximate 
form (North, 1975a; North and Coakley, 1979) 

S(x) S 1 + $2 P 2 (x) , (3) 

with S, “ -0.477, and ? 2 (x) is the second Legendre polynomial, 

P 2 (x) ■“ l/2(3x^”l). With this approximation, S(x) is a parabola In 
X, having zero derivative at the equator (x “ 0) and falling to a 



value of 0*523 &t the pol® fe ” We t.ote in i>asslag that a; is a 
convenieat variable to ase ia sooai average applications because dx 
is proportional to the area of a latitude strip and therefore the 
area aver?4ge of q(x) for th® region apannad by da is given by 

<i “ q(s) dxfix « ( 4 ) 

Using a value of 0.30 for a (Ellis et 1978), we eoapute 

254.6 K for the Earth’s radiative teaperature. Clearly, this 
is aiuch colder than the observed saa level average tenperature 
which (for the whole globe) is 287.4 K. The major part of this dif- 
ference is, of course, due to the so-called greenhouse effect of the 
at»sphere to be discussed Isxer* 

Let us now cosspute the fundaMntal sensitivity paraiaete^r, $ 
defined by 

0 dX 

o 100 do ^ 

o 

where the subscript o on T refers to the global average vaL 4 tie « 0 is 

o 

a laeasure of the change in ^obal average temperature due t-to a 1% 

change in the solar constant. For all climate models 0 isv the first 

o 

quantity to cot^ute because the sensitivity of the model toy any 

perturbation is rotxghly prcfOTticnal to 0 (cf. Sections 7 rand 8). 

o 

Vith constant the sinp3.c raodel defined by (1) give.s 


0 (black radiator) T^/400 ® 0®63 K 

O R 


( 6 ) 


Tills number repx^esCTts the sensitivit:; of a system with no ff^edbacks^ 
and it is a ©taiadard for with all cXlmte ridels... 







Tlje sensitivity of the actual climate is influenced by a oiriad 
of feedback processes (Schneider and Dickinson, 1974), For exaiEple, 
changes in planetary ten^scrature night change the emissivlty or reflec- 
tivity and thereby modify the sensitivity. Many of these agents and 
linkages have yet to be Identlfiad. Of the known and suspected feed- 
backs tsany of the details remain a oystery. Some of these uncertainties 
are discussed in Section 4, Here we introduce some eliisple models to 
illustrate the potential impact of known feedbacks on climate. 

It Is well to note here that the value taken for the solar constant 
a . 1340 W m”^, is different from that used in Section 4 (1360 W m ) in 
order that we not have to repeat published calculations. In fact, various 
modelers have used different values for this parameter in their sensitivity 
studies. Generally the sensitivity results are not sensitive to the 
present value of the solar constant. Errors Introduced in this way tend 
to be partially compensated for in adjusting other unknown parameters to 
force the model's unperturbed climate to fit the present climate, a 
process sometimes referred to as "tuning.'* Incidentally, the exact value 
of the solar constant and its constancy In time are the subject of con- 
siderable experimental activity both in rocket measurements (Willson 
et al., 1980) and satellite measurements (Hickey et al . , 1980), Tlie best 
current val e from the Nimbus 7 satellite is 1376 W m , 

a. Global models with feedback 

Feedbacks affect the sensitivity through their Influence on the radi- 
ative fluxes absorbed and emitted. Often the net effect oi the feedbacks 
is inferred from empirical data. We begin with the flux of infrared 
radiation emitted by the earth. 

Budyko (1969) suggested that the infrared radiation to space can be 
represented as a linear function of the surface temperature T (in *C) 

I - A + BT , 

where A and B are constants deduced from observations. Based on data 
from the northern hemisphere (North and Coakley, 1979) we find that 
A » 203.3 W tmd d “ 2.09 W gives the best fit between the 

fiujces caiculafed using (7) and those observed. The energy balance 

©ay be written " 



I 



A + DT 


Q(1 “ «p) 


( 8 ) 


* 


^je^e Q Is or /4. Using (1 ~ o ) - 0.70, we arrive at T « i4,97®C, which 

O p O 

agrees with the northern hemisphere value (14,9®C)o The coefficients 
A and B take into account average cloudiness conditions, the effects 
of infrared absorbing gases and the variability of water vapor» For 
cossparison a linear escpani&ion of c?(2?3 + T)^ would lead to "black 
radiator" coefficients « 314*9 W and “ 4.61 W a 

For constant albedo, the sensitivity of this "greenhouse model" 
is 

A + BT 

3^ (greenhouse) * ** 1.12®C • (9) 


We deduce that the presence of an atiaosphere increases t^e sensitivity 
of the climate. This effect is referred to as a "positive feedback" 
since it increases the sensitivity over that of a black body 
radiator. 

Let us examine the reasons for this positive feedback. Consider 

a planet surrounded by a shield at temperature T^. We will assume 

that this "atmosphere" does not absorb solar radiation, but perfectly 

absorbs infrared radiation. We imagine the shield to be in equilibriiim 

so that it radiates (net up and down) at a rate equal to the rate at 

which It absorbs. As a result, T, is related to T by 

1 o 

ax ^ “ 20T: ^ 

o I 

Similarly, if the surface is In radiative rquillbrlum, then is 
given by 



10 


. m, 

nom <10) w« deduce that i» gysater than by a factor of 2^^^ 

(8^ 1.19). rurth«n®oire. hy cosibining <i0) and (11), m sm that X 
la the effective radiating tee^erature of the planet* For « 234 K 
wa compute 7^ « 302 K. Ihua a single black ahiald approxiaatea the 
greenhouse correction of the earth's attuosphere. Co«&lnlng (10) and 
(11) and diffarentiating vlth respect to Q, keeping tha albedo 
constant g va obtain 

Bjj (black ohleld) * T^/400 » 0.76 K . (12) 

Although the black shield approxioatfcs the attaosphere'a greenhouse 
effect, it foil/* to account fo? the increaotd sensitivity when the 
etSMJsphere ia present. 

The Boat probable reason for the enhanceaent is the variable 
concentration of water vapor. Water vapor Is a significant absorber 
in the infrared. On the average, the amount of water vapor in the 
atmosphere increases as the teaperacure increases vllanabe and 
Wetherald, ^967), To allow for this increase, we should allow the 
maaber of black shields to Increase with terafserature. The empirical 
coefficients in Budyko’s formula prcsusMibly take this effect into 
account. There are, of course, other feedbacks that affect eslsslon 
and their influence is also reflected in the coefficients. In the 
next section v« shall return to the Budyko formula and discuss its 


validity* 



Left \m ®cw «satain« the Sce**cep albedo feedback* If the Earth 

had 00 ic» or srom but: it had 50% cloudiaetts* the local co-albes!© 

(I * albedo) %?o«ld ba spjsroalsately a^(s) » 0.70{ »is*ll«rly» on 

tc««cov 0 >ad plsast would have cO“«lbedo a^(s) « 0.38. Using (S) we 

coi^ut® V ■ 15.0®C for the ice~fre© planet and •» ~36.4"C for the 
o o 

ic«8^cove plmete Of couira# cha aurth lim betvean aKtrema* 

hist^ m pTOpoae a end® i»d@l for the ica cap aisa* Stippoae 

that the edge of the ice cap is denoted by x • x^. T!\cn suppose 
Kg « 1 for > 15*. »g • 0 for < -15*, and In between 

Kg " 1 (T^ - l$)/30* We Bay then use for 0 ** x < Xg. and 

for X < X < 1, From (3) the planetary co>albedo would then be 

o 

given by 

“ " 1 - Op 

■ »1* [». + I Sj(., - iiolx,(y). (i3) 

The notation for the planetary co-albedo will prove useful later. 

H fx <T )1 ia represented by the solid line in Figure 1. To find 
solutions for (8) we plot the outgoing IR divided by Q, represented 
by the dashed curve. in Figure I. For the present solar constant, we 
obtain three roots. The root labeled X corresponds to the present 
cllaate, Root II is an intermediate climate with the planet having 
about 30% of its area covered by l,.e, and Root III Is on ice-covered 
planet. Such lajltlple solutions of sero-diraenslonal models were 
noted by Sellers? (l§7&), Crafoord and Kallen (1978) and Fraedrich 
(1978), 


Xn Hgure 1 m have a qualiteelva picture of th« es»d«l ciiaataa. 
Xf Q is 4«creasftd» the dashed curve is scaled t^ard en<;d louts XX and 
XXX Gisrge cose the r and then disappear* leaving only the deep freere* 
loot Zn* if the solar constant is increased* loots XI and XIX 
coalesce and ve are left only with an ice>ire« larth* l^oot X. 
rigure 2 shows aa a function of Q. 

The esas^lc just given illustrates the rich atruccure of even 
the Ice feedback Kiodel* The eulti-ple branch structure shown 

in Figure 2 persists in laodels with latitude dependence end even 
exists In CCMs (Uetherald and Manabe* 1975). 

b. One-dlmenslonal taodela and transport 

After the cero-dlcwnslonal oodels, the next laodelis to be studied 
include latitude dependence. TI>1 h additional degree o-f freedoo 
forces UB to consider the horizontal transport of heat by the geo- 
physical fluid system. In the treatment of this transport we will 
have to sake drastic Idealizations to keep the aatheiaacics manageable. 

In the latitude-dependent ludels we assume that the rate at 
which heat enters each infinitesimal latitude belt during the year la 
exactly balanced by the loss rate. The individual ceriaa considered 
are schc^mntically represented for the istrip: 


(net hotiKontal tranf^port out). ^ (infrar^ ■ 

^ " (14) 

■w (solar * 


A coimon fa^ctor* the area of the atrip^ tov he cancelled throughout 
to that t\w regaining termr. have tmitet energy per t^nlt area 
per imlt ti» 0 u 






M ••«>« C» .U, U'i,^ ,1.0. 



n 


proceed vith the crucial ei«^llfylng ««8ffiu^ticms that each 
tts'tm in (14) isay ha represented as a iuactioa of the aenally av^tregad 
sea I@% , i (1000 ;d>) tessperature field* T(»j^). Once the pssrsassterisatlon 
fofffloles for the individual teres in (14) ora found* we ®iolva th® 
ayoteia of equations reprasented hy (14) for all indices 1. aiml** 
taneoualy. It is usually nf^casasry to is^osa a boundary jconditlon at 
tha poles* since tsathettatlcal eystasss like (14) often have an 
unphysical solution (Irregular solution) that diverges at: the poles* 

By Itaposing the condition that the flux of horlsontal he^t into 
the pole vanishes* we systematically eliminate such spurious model 
climates. 

To aee hoy the model behaves* let us exat&ine a few e^^.'treue 
cases. We adopt the Budyko formula (7) for infrared radMatlon at 
each latitude and we adopt an ice cap parai^terizatlon tfeat is also 
due to Budyko (1969)* namely that the ice cap edge extendls to the 
c»an annual laotherm 

T(x ) ■ T «■ -10*C . (15) 

S 8 

For KiTsipllcity , we take the co^albedo 8 ( 25 ^ to be dl£?<i;^otttin\iou0 at 
the ice cap edge* Note that mxtM come sotssehow from Ql5} rather 
than from tha llne"'r relation with that was adopted im the flobal 
average »del* 

For the that the horizontal transport Ijs the 

planetary surface be isothemil* This Mfdel ta the 

glolniil average? model that the betveen *J5®C IS* ®hown 




14 


In figure X ia coaspressed eo a vertical line centered at ••lO'C. The 
solution branches are easily eho«m to be given by the dashed lines in 
figure 2* Tfie physical interpretation la obtained by lolloping a 
quasi-static change in Q de?m Branch I of the solution. When the 
dark (Isothsmal) planet’s tetaperature reaches -10*C* the planet 
suddenly turns white and its new equilibrium tei^ereture isust be 
about -34*C. 

For the case with no transport the enex^, Helance equation 
becomes . 

A + BT(«) - Q S(x) a(x, »^) . (16) 

at each latitude. We assume further that at x * x . the co-albedo 

8 

a(xi X ) bacomes the average of a. and a.. 

B il £ 

Applying (16) at x - and using (15) we obtain 

s 

Q(Xg) - (A -t BTg)/S(s^)a . (17) 

with a “ (a^ + a^)/2. Expression (17) gives the solar constant Q 

needed to maintain the Ice line at a particular latitude given by 

X . The curved line in Figure 3 is computed using (17), Ttie lower 

flat line (x » 0) is obtained by starting with coasplete ice 
8 

cover (Q/Q^ 0.8) and raising the ©olar constant quasl-statically 

until the equatorial tciciiperature reaches which occurs at 

Q/Q^ ** lel7)a Similarly the upper flat portion la obtained by 
lowering the Kolar constant from Q/Q^ 2*0 on the lce--free planet 


until the pole xreaehes -10^C» The peculiar solution lines emanating 




fross £ha curve in Figure 3 are cabled cvi-8s> by ®ath«?wtlcian9. Ve 
shall return to the structure of these solutions later* 

The sx>dels discussed thus far present a paradox. If the sun's 
luninosity has risen fresa a value of 20 to 40Z lower than its current 
level, as virtually all solar evolution theories indicate (llcuctan and 
Sood, 1977), Chen why isn't the earth covered with ice? The taodels 
indicate that if the solar constant is raised qus8i~static:.lly from, 
say, three-quarters of its present value up to its present value, the 
laodel stays on the lower branch of the solution cun'* whether the 
transport is Infinite or zero, Furthenaore, we shall see th.at tlila 
Teeult holds in every aodel studied in this paper* We are not in a 
position to resolve this solar evolution-climate model contradiction, 
but conjectures have been ventured: not enough nolsrure is available 

at the right places to generate total ice cover; other negative 
feedbacks are present, perhaps related to cloudiness change (our 
calculations are based on present cloud cover); and the cot^»osltlon 
of the atmosphere was different in the past such that a larger green- 
house effect prevented the ice cover (Sagan and Mullen, 1972; Owen 
1979; Hart, 1978; Budyko, 1977). 

Clearly the no-transport model bears little resemblance to the 

earth since a solar constant 70% greater than present Is required to 

push the ice cap b/ick to » 0.95 (Its present Iccation), With 

Q 1,7 times its present value and • 0.93, the planeta^Ty average 

teiape rfsture b«coaisp^s ^ 93^C* Asidf? from the unrealistic values 

obtained for T end x , (16) suggests that T(x) has a discontinuity 
0 8 
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of She order of SO'C ®t » " Transport, of conrae, tends to 

6f^oth this discontinuity* 

Before proceeding with a finite transport Bodel, let ws look at 
the tonal average tesiporstures obtained for the ©stress cases of Infinite 
and aaro transport. For the absorption, we take observed (northern 
heaisphere) t»an annual values which are reasonably well-represented by 

a(x) -80 + 


with a * 0.681, a» * -0.202 (North and Coakley, 1979). The planetary 
o 2 

average tetMperature is obtained by integrating (16) over »; 

A •+■ BT “ Q f dx S(x) a(x) 

° Jo 

Using the constants specified earlier we obtain “ 14.97“C, Figure 4 
shows T(x) (solid line) for the no-transpert, no-lcc-feedback case 
computed directly from (16). For comparison, the figure also shows 
the Infinite transport case (dashed line) which yields an isothermal 
14.97“C planet, and the teiaperature observed for the earth. The 
transport has the obvious effect of warming the pole and cooling the 
equator. It snalntains, however, the roughly parabolic shape of the 
no-transport case* 

Consider the addition of a transport term to the heat balance 
equation. Tlie geophysical fluids transport heat through their ciaan 
and transitory (eddy) motions. If we average the attaospherlc 
velocity field through, say, a month and around a latitude circle, we 
would sample wny statistically uncorrela ted eddy processes. A 


evsw. 




first mdel for the atssoapher® is a geophysical fluid having a random 
velocity field* U@ edght taka the heat content es a passive or .»lar 
being carried by the fluid* B(}ual atsotsnta uf fluid are directed 
north and south across a latitude circle in any Interval, On the 
average, heat is carried from vam areas to cool by an amount propor- 
tional to the gradient of the temperature: 

heat float « 

where C is the heat capacit> per unit area. The proportfconailty 
coefficient in this simple model would be a diffusion coefficient. 

The amount of heat per unit time per unit area leavlmg a strip 
is the divergence of the flux which is proportional to or in 

our notation 


__d 

dx 


D(1 - x^) 


dx ' 


( 21 ) 


where o may be a function of x and must be thought of as aa free parameter 
to be adjusted effiplrlcally. 

Obviously (21) represents a gross overslmpllflcatlom of the trans- 
port process. Ihe most evident omission is the mean circulation in both 
atmosphere and oceans. We shall consider more general mo«5els after 
first studying the model defined by (21). One advantage ajf the form 
(21) is that it coi. tesponds to a physical analogue, namely,!;, heat con- 
duction and, therefore, physically realizable solutions arse guaranteed. 
With transport the energy balance equation becomes 

” dx D(i - s^) A -f BT(x) » Q S(x) a(x, x^)) . (22) 





Tha bomidary ccmditicms for ©ywsatric solutions ar® by 

- D (1 « dlCsI 

dx 

Note that the eyateia is nonlinear du© to th© ica-tilbedo foedbaclc. 
Th© nonlinearity leads to the multiple solutioi^s that ^ have already 
encountered with the aifspler siodels* 

Before solving the nonlinear system^ let m insert observed values 
of the co-albedo a(x) in (22) and see how the coi^puted temperature 
field compares with those shown in Figure 4* First talu D to be a 
constant independent of %• Since the Legendre polynomiidLs are eigen- 
functions of the diffusion operator, 

d 2 

- ^ (1 “ X ) — « n(n + 1) P^(x) , (24) 

and since individually thsy satisfy the boundary condlttcis 

(23), the even-numbered P^(x) form a convenient basis a©i for expansions# 

Consequently, we expand the surface temperature, 

*<’■> • S • «5) 

n 

even 

and insert the expansion into the energy balance equation (22) ; 
noiltiply by Pj^(x) and integrate over all values of x frosiO to 1. 

Making use of the orthogonality relation 
fl 6 

dx P^(x) PJx) « ra,n even , ( 2 fi) 

we have 


- 0 . (23) 

» - 0,1 
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LT+6 A»QH , 
n tt on ^ n ' 


( 27 ) 


where 


L » n(n 1) D 4* B , 
n 


( 28 ) 


and 


H (x ) - (2n 1) 

n 8 


P (x) S(x) a(x, X ) dx . 
j n 8 


( 29 ) 


Even though we Intend to use the observed a(x), we hove retained x^ 

In (29) for later applications. Since the a(x) are given by (18), the 

H are determined, and as a result, (27) is easily solved: 
n 


T - Q H /L - 6 A/B , 
n ^ n n on 


(30) 


For n « 0 » 


Tq - (QH^-A)/B . (31) 

which is equivalent to (19) • As before, is the planetary average 
temperature a The form of (31) holds not only for the diffusive 
transport (21) but also for any term that is the divergence of a 
flu 2 C which itself satisfies a zero condition at the end points 
(equator and pole). Using the coefficients specified earlier, we 
obtain as before, 14.97*^0. 

We turn now to n « 2. From (2.3) and (2.19) we compute 
« “”0.500. In order to compute T 2 , we must know D, which was 
left as an adjustable para«ter. By taking the observed value for 
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(-28.0*0, we confute D • 0.649 H eT^(*C)’'^, or in dissnsionless 

fora L«/B “ 0«310« Thie approsdlMtion is called the two-iaoda or 
parabolic approscimation to T(x) (Isagstrom^ 1928) t 

T(k) « *^0 

It is plotted as the solid line in Figure 5 along with the observations 
Of course » the agreement shoro in the figure does not confirm the 
model but rather illustrates how well a parabola fits the data, 

I^et us now consider n « 4, From observations (Ellis ^ al, ^ 

1978) we obtain »« 0,022, and using the previously cos^puted value 
of D, ve coispute « 0,5®C, The observed value is « «3.5*C# 

The model gives a value that is almost an order of magnitu<^e too 
small and has the wrong sign. In fact, the higher mode amp^litudes in 
the model are so small that the two-mode solution is close to the 
exact solution (North, 1975b), Perhaps the most prominent feature of 
the observations that is not reproduced by the model is the^ flattening 
of T(x) Into a nearly Isothermal band in the tropics. This flattening 
is probably due to mean taotlons in the tropical atmosphere (Hadley 
cell). The mean szkotions are more efficient than random motions at 
suppressing temperature deviations, Lindzen and Farrell (1;977) have 
discussed methods of allowing for the tropical mean motion©;, Tneir 
motivation was based upon results derived from a more compir}:ehensive 
model (Schneider and Lindzen, 1976) in which it was found tihat near 
the equator T(o) - T(x) « In terms of the Legendre cxp/anslon, 

temperatures near the equator require contributions f rom (x) • 



IB an ^ hoc way we caa Issitete the efficiency of theimal conduction 
in the tropicB hy «no»ins the diffusion coefficient to depend on r; 
e.g.» IK») " *•* patssaters, 

end D^. that can h® adjusted to ^ve the observed values of and 
T^. In agceemnt with Undren and Farrell, we find D(s) large near 
the equator and smII near the pole (G. R. Korth, P. B. dates and 
R, F. Cshalan, unpublished). Although auch tricks suffice to correct 
the value of T^, other processes, such as varying relative humidity, 
lapse rates and cloudiness, that would affect the Infrared caJssion 
could also contribute to T^. Further discussion of the lifflitatlons of 
th® diffusive epprosiaation are postponed to Section 4. 

Let us Interpret the expansion (25). Each succeeding term in 
(25) contains inforaation pertaining to smaller and smaller spatial 
scales. The first few terms give us the gross features of the 
planetary climate. is the planetary average; (3/2) is a rough 
measure of the pole-to«*equator temperature difference; higher order 
terms reveal features at finer spatial scales. Hence, the spectral 
method of solving clltsate models provides a fraj^avork in which 
aathemstical (or numerical) technique goes hand In hand with the 
concept of a model hierarchy discussed earlier. 

The two-mode approximation just discussed was derived for dif- 
fusive heat transport with constant coefficient. The form of the 
transport, however, need not be restricted if (25) is truncated at 
a » 2. For Budyko <1969) instead of using the diffusive 


tex*m (21) 
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YETW - T^l 


(33) 


^ere Y is aa rapirical coefficient siaillar to D® Budyko^s model is 
thiss given by (22) with (33) replacing the fi rst term* In the 
two^mode approximation the mod«:is are identical. By aubstitntlng 
(32) into (21) we obtain 


« 6D[T(x) - T^] , 

v;hich is identical to (33) with Y ** 60. While the two models are 
Identical in the two-tcode approximation, they differ in higher 
modes. 

Budyko’s iiaodel can be solved analytically (Chylek and Coaklcy, 
1974) by roughly the same method we used for the no-transport model. 
After some straightforward algebra, we obtain 


A + (B + Y) Tg + ir A/B 

" S(x ) a + Y H (x )/B 
s s o s 


(34) 


Note that if Y “ 0, (34) becomes the solution for the no-transport 

model (17). The solution curve for this version of the Budyko model 

is shown In Figure 6. The twigs at Q/Q « 0.99 and 1.13 come about 

o 

in the same way as in the no-trsnsport model. Figure 3« 

Although there is a clo&i^d form solution to Budyku*ij ^del, (34), it 
is instructive to investigate the spectral properties of the model. 
Expanding the temperature field as in (25), and inserting the eitpanslon 
into the energy balance equation, we obtain again (27), but with 
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given by 

- Y + B - 5^^ Y • 

I„ the d«£uh<v. «.1.1 •£„ a l/t,. £<.». It D > 0, th. hith-t «»l.e 

(„ > 2) «M MtohSly suppressed. lu the Buiyko ■»d«l. su thei other 
hand, the higher nodes are not suppressed. As a result a dls»c 
tlnulty in .(a. a.) lead* to a discontinuous t«®.ratur. llel.d 
In the Budyko model but not In the diffusive model. 

Ihe time dependence ot these models also differs. If »e add 
heat storage, C3l/8t to th. energy balance equation and eapamd th. 
tmeperature fields, then ». obtain (27) «lth the additional norm 
Ci„ on th. left-hsnd side. If the soUr constant 1. suddenle. sultched 
off, the Indlvlducl node amplitudes decoy exponentially vlth tlnm 
constant C/I . In diffusive models features or anomalies that hove 
small spetlel scale, decay rapidly, vheress In Budyko-type models 
small space scale features decay at the rate of the planetar-y scele 

features (n » 2)* 

ut US no. consider sosm generalisations of the diffusive model 
in th. t^emode approrlnatlon. Suppose D depends on latltudie, 

B(«) " nwW. Ih. n - 0 relation. (21) 1. unaltered and th» n - 2 
version of (30) Is also unaltered provided ve replace » In <28) 


by D’» which for n “ 2 is given by 


n 


^ dP,(x) ' 

- -5 nj 

*1 ^ 
o 

(1 - 


Since the Integt.l In (35) Is a constent, to be computed once end 
for all. It may be absorbed Into the phenomenological constant n. 
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Th® constens >1 is to b« «8 before* ao Sthat in eh« Kwo*S8od® 

epproKimation D(k) » nw(s) i® feimlly equivalent coa®t®at D 
(Uorch, 1975b). 

Thl© formal tquivaleaee is also retained for ©tbsr trenaport 
®odel«» For esaisple* S«llerB (1969) iaclwded & msiiaasi circuletloa terra 
given by 


vcsl a - (36) 

vher® V(r) 1 b an datcraissad m^v©c£iv® %wIocity fields M 

«bove» in the appt© 2 dLmti^a#this t®ra my? be absovbod into 

the coMtant D® A^iother form of interest Is K?^T(fe)» which could occur 
as an added tora in a turbulent (rsndois) atmospheres^ ^del (Norths 
1976; Kails, 1976)# Because ^P^(x) «» tn(n 1)1^ " P^(k), the fomal 
equivalence to constant diffui.ilon in th® two-®y>de ^^i?.ppro^mation8 i® 
maintained^ 

Tlie iomaX equivalence of so many different in the two- 

i25iode approKimation explaina why various published saaodels give the 

saw qualit&tlve results® Of course, it is xisicei&&>.,i^ity that the two 

mode%'5 represent a suitable fit to observations* lit i& &1 bo necessary 

that the two-TOde nolution approxiMte th® exact snf>Xution when the 

nonlinear lc«-slbed feedback Is included* To obt^^ jiu an exact 

solution, we my solve the system (22) analytirallvy* Here we retain 

thn % depiendence In H (x ) (27)* Dividing (27) thv?rough by L , 
n ns n 

rauitiplyJEng ty P (x ), Bussing ov«r n and using (1^5), ve obtain 

n S 


V 

A/St 4- T « Q / 
tx 

even 






n 


or 


A ST 




(37) 


even 


The®® solutions tmve bean given fe? end Sueres mU) end by 
Roreh (1975a,b)^ 

Rote etmt the fomal «ol«tioa ( 37 ) also holds fm the Sudyko 
^del or even « KV y mdel If the ere apprspri«t«% ss^difled 
{since It .8 fonaaUy correct for any a(x.«^)l, the solution is not 
restricted to e discontinuous elbedo. He could use a aw,ethed albedo 
at the lc«“cap edge such as 


«(st.Xg) - ®o + Cj tanh {(» - *^)/wJ , 

as will be discussed later. 


(38) 


M 


of tri^nsport 

To generalise ( 23 ), let the energy balance equattoi be given 


by 


bfT] ■¥ A(x} « Q S(k) a(x,x ) 


(39) 


where I. la a linear operate- (Cahalaa and North, 1979).. Tne te>del 1b 
cosaplcted by the ice-line condition ( 15 ) and boundary cmditlons 
similar to ( 23 ), For constant diffusion 


L » “» 4" (1 “ ^ B . 


( 40 ) 


In general, however, L might include other Integral or -rnfferentlnl 
opetatot's, ■ 










Ua Bfesll ttsa f»sacti®a 0 (js,x*)) for 1. axlstSt 

&Q that 

t^G^Csc»s*) - 

wh®r« 6(s - s’) im th* Mr«c dalta^wctlo® ss^ tha &i&mattpt s cm L 
indicataa sls«t ossly th« s varla&le 2» effected &y L. 'E«® eoiutlon, 
to tisa linear iiirasogeaaoa* fsss&lej®. 


“ P(X) 


( 42 ) 


where 0(») is a kaowa iwctic® Is saiqualy givsa by 
J(x) - G^(x.k') p (x*) da* 


( 43 ) 


provided L ead the «iss4xl4ited h-minfery conditions lead a tmique 
Inverse. {For properties of L n«^«sary for the exlstenvce and- 
uniqueness of see> for feurant and Eilbcrt* 1953) ♦ The 

Green*® function for cc»tent diffmion in given by 




(2n I) P ix} P (x*) 
, ?> n 


eve^s 


(44) 


(A closed-fora espressics for (44) I* derived ia Appendix 11 ,) 
For Istltude— depeisdent diffusim 









we us« the « i gen functlos:^^ deifined by 


( 45 ) 
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4 5 d f (a) ] 

- S ^ J “ Pjjn) , (46) 

wItSM fchs era ®lg«v®itjss8. TS'sa syste® (46) Sotm a St«r@~Llouvine 
® 78 £®a (Cmtmt md IliSfeert, 1953). prwided D(sj) > Q mu U w«ll- 
bahavad, Th^ g^(a) propertioMl to l»^<s) for » » constant. 

The aS. 0 o are orthogonal and can be aorsaslised ewch that 

1 ^ • (L7\ 

® im v47) 

Further^re, the lowest eigenvalue U - 0 and f (a) « 1 and for n 

2 o » 

Wjj ® n . The Green's function for the ayatesa la given by 


C'”’ <»■*•) ■ X f„(») *„(-•)/»„ . 


wheris B# 

The Green's function G^(x.x^), may be Interpreted as th* thermal 
response of the model to a localized heat source at x •• x . Figure 7 
shows a graph of this function for the constant D diffusion model. 

For large D the heat Is smeared outj for D tending to aero the response 
tends to a spike at x “ x^. One ii^iortant feature, to be used later, 
la that G^(x,x^) is positive definite. 

The formal solution of any model possessing a Green's function 
Is straightforward. Using (42) and (43), we convert (39) to the 
Integral equation 


T(x) - / G^(x,x») {0 S(x') a(x‘,K ) - A(x’)]dx’ , 


As before, we satisfy the Ice line condition (15) by setting 
In (49). Solving for Q we obtain 
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QCsJ 


f 


Tg *»- J 6^(s,»s') A(jt*) 


r 6 S<s*> aCs'.Sg) d»' 

i ek 


(50) 


Tt^la solution reduc«» to tb« special cases (34) «asi (37) **00 the 
appsroptlste substitutions ate isade* 

Figure 8 shows the Q(s^) versus curve for the constant D 
diffusion Eodel* computed using (37) with teroa through r« " 6* The 
exact solution (North, 1975a) differs little fro® that shown in 
Figure S. The most striking difference between the diffusive model 
and the Budyko model is that the twigs in the Budyk® snodel split 
open Into two branches in the diffusive model. Just below x^ ■ I, 
the slope changes sign, and dx^/dQ tends to aero at the cusp where 
the two branches meet (Dratln and Griffel, .1977; CMialan and North, 
1979), If Q is Increased quasl-statlcally from Q^, Figure 8 indi- 
cates that the ice cap first shrinks to the point ^5ei‘e the slope 
changes sign, and beyond that presumably a Jiuap to m ice-free 
earth occurs. Let us examine whether this sign change and the 
associated cusp ere physically realistic. Lin (19.vV) -and North 
(1975b) noted that because series (25) converges ra?ldly, suodels with 
nonlinear diffusion coefficient (D « dT/dx) are easily solved. The 
two-mode results for the nonlinear model arc virtually identical with 
the linear model but retention of the higher BK)dcs 3a the nonlinear 
model removes the cusp at x^ ® 1 shown in Figure 8* A picture more 
like Figure 6 emerges. 

n« cusp can also be removed by smoothing the albedo at the Ice 
lln®. Cahalan and North (1979) have expericenteJ ii^th the aavoothed 


c’, . i:- ^-0- 
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£osia (33) choosing Cj and Cj to wtch tha valutas off the corresponding 

step albedo but allotdng w, & s^aauta off tha smoothing width, to 

diffet ffrc® sexo. As the width of tha cmoothlng increases, tha cusp 

disappears sod eventually the slope bacoasss positive near a •» 1, 

s 

Coakley (1979) also laentions this affect* For smooth albedo, the 
twig in the Eudyko model at a » 1,0 is also ressoved. Cahaian and 

cl 

Hor&h shorn that swh feahavior follows goaatally firom (49) and 
(50)* We m^t conclude that the appearance of the twig or casp ne&r 
» 1 in the ©issple mdela is probably ^physical and ts^rely an 
artifact of the Mthematicaily convenient but physically imreallstic 
step-function albedo^ 
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* Se®©otasl TOdeX® 

Xn estendiiig th« one-aiaenslonaT. cUeate ®od«l to Sacl«5* 
seasonal cyclss. we wish to ®t«<Sy the sensitivity of the eodel to see 
1£ it differs from that of the corresponding «ean snnual ®odel. In - j 

other words, do ceasonal changes effect the nean annual cUmte? 

There have been several attea^ts at developing seasonal ©odels, 
notably the energy balance aodels developed by Aden (1962) and 1 

Sellers (1973) and the general circulation nodels developed by Manabe 
and his co-workera (Wetherald and Manabe. 1972; Manabe etal.. 1979). j 

More recently, there have been studies by Thotapson and Schneider > 

(1979) who seasonallzed the Gal-Chen/ Schneider version of the Selle , 

«,del (Schneider and Gal-Chen. 1973; Gal-Chen and Schneider. 1976); 

Ramanathan et al. (1979) who seasonallzed the Budyko model; and North 
and Conkley (1978. 1979) who seasonallzed the simple diffusive mode. 

(North. 1975a,b). The model studies so far have similar qualitative 
conclusions so that it will suffice here to study the diffusive model 
as fin exataple. 

As before, the method of solution (Legendre polynomial expan- 
sion) , is intimately connected with the philosophy of the model 
hierarchy approach. The development of the diffusive seasonal model 
is similar in spirit to the approach taken some time ago by Fritz 
(1960). Our em, -basis here will be on using the model to unde 
th. «.ln £.awt.s of tt,. .o..on.l cy=l«. *s f«t .. possIbU, vc »IU 

use analytical tiiethods. 

Mo bogltl by oitoiiiliilns ee»"n»f , 

foporotur. t(y.t), >1.0 infroboa flux tho co-olbodo o(x,0. \ > 
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tisa hsat per mit arsa raachins th« top of th« at8»js?ph®x® QS(s»t)» 
north an4 CoeUloy (1978, 1979) have ehcnm that a convsaf sat rep.*®- 
8®tttatloa of the fielda is to first ayraetrire thea so that data froa 
only one hemisphere is used; the other hemisphere is given the saae 
data, but lagged by 6 months. After syra^trisstion, this data are 
fitted to ®t®pl© fosmlae 

P(*,t) " cos 2Tft + BjjSin 2irt)Pj(x) + $‘ 2 P 2 («) » ^^1) 

where F(x,t) can be any of the four fields. Each field Is then 
characterized by the four coefficients in (51). Note t&st the oean 
aiuiual fields are characterized by F and F, while the ais^lltude and 
phase of the seasonal cycle ®ay be cosaputed froa end ^^.i* 
for the four fields using the ByMaetrIxed northern hemisphere (SN) • 
syMoetrixed southern hemisphere (SS) and global data haws been per^ 
for^d by North and Coakley (1979); results for SN are l:isted in 
Table 1* Figures 9-12 show exas^les of the curves for warious 
seasons along with observations. From the curves ve see that (51) 
captures the gross features of the fields, but misses the fine 
structure which is particularly noticeable at low latit^^^des. In 
keeping with the approach taken so far in this review, we would not 
expect simple models to work on scales beyond p2(3t) in C51)# Indeed, 
attempts by No^ith and Coakley to go beyond this failed. 

As in the mean annual wdels, we adopt the Budyko form for the 
emitted infrared radiative flux, 

I(x,t) A 4- BI(K,t) 


(52) 



mth A • 203.3 and B « 2.09 W (Moreh rad Coakley, 1979), (52) 

gives the best i±t to the TOde ©j^litudes listed in Table 1 for the 
emitted f I«es and temperatures of the northern hemisphere. The fit 
to observations is ehora in Figure 10. 

We construct the seasonal model along the lines of previous 
sections. In the seasonal mdel, however, a storage tent, C(k,<^) 3T/3t, 
tausc be added; C(x,<J)) is a latitude and longitude dependent heat 
capacity per unit area. The forcing is given by the four amplitudes 
in the ©ode ejcpansion of S(x,t). The seasonal response is repre- 
sented by the four amplitudes in the ©ode cKpanslon of T(x,t). A 
linear climate model would connect these two fields with a four-by-four 
response matrix. Linearity is suggested by the absence of higher 
harmonics in both the forcing and response field; (North and Coakley, 
1979). The most general linear model with four components consists 
of a response matrix with 16 Independent components. We shall see 
that most of these elements toay be taken to be zero. 

With the storage term added the energy balance model becomes 

C(x,(|i) 1^ - DV^T A + BT - QS(x,t) a(x,t) . (53) 

at 

As an Idealization of the northern healsphcre, we take a single con- 
tinent with coastlines running along taerldians and with an area 
spanning 40% of each latitude belt. We separate tne continent from 
the oceans because C(x,({>) is only about (0,16 B) years (C^, ) over land 
and about (4.7 B) years (C^) over the ocean mixed layer (75 o). 

If (53) Is integrated around a latitude belt, we obtain for 
the lasu! mass (after dividing through by f^^, Che fraction of land area) 


left 


3T^ <x»t) 
dt 


^ 3 2v a „ "o 8T 


right 


( 54 ) 


A + BTj^ •• QS(«,t) aj^(s,t) , 


where T, is defined s» the average teaf>erature over land i» a latitude 
L 

belt. The gradient tera Is the difference between land and water 
tesnperatura divided by an effective angular distance over which the 
tesiperature change effectively occurs. The whole tero stay be inritten 
as 


^ (Tj, - V . V > 0 . <55) 

la 

vhere V is a new adjustable paraiaeter that accounts for the land-sea 
interaction. 

An equation analogous to (54) can be derived for T^(x»t), the 
temperature over oceans* The ronal average field is the weighted 
average given by 

T(x,t) - fj^Tj^ + f^T^ . (56) 

This kind of dccosaposltion was first used by Sellers (1973) • 

If w«* if,abstltute the truncated Fourler-Legendre series (51) » for 
T, and T into the resulting equations, we obtain 

L W 

A + BT - QH , (57) 

o o 

C — + (2D + B) T,^*” + (T T - QH, (58) 

L,w dt 1 fL^„ 1 11 

(6D + B) T, .« QH, 


atid 
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where 


H^Ct) 


H^(t) 


H2(t) 


~ { 23 q + a^<t) Sj^(t) + 02826^22 + .« } . ( 60 ) 


2 ®101 **■ ®1I0 %i 


+ ®l<t) S2G222 + ••< 


) . 


( 61 ) 


T ^®2®220 ®0®2®202 '’’ ®2®2S22**” ®i(*=)Sj('fc)®21l‘^ ***^ • 

( 62 ) 


and 


'ijk 


-1 


Pj(3£) P^(x) dx . 


(63) 


The relevant coupling coefficients are « 2/3, C^22 “ 

^121 * S 22 " (60)-(62) Sq - 1; a^(t) and S^(t) 

represent the first harmonic contributions to the coefftiiiBnts of Pj^M 
la the expansions of a(x,t) and S(x,t), For simplicity ^ have assumed 
the co-albedos for land and ocean areas to be the same* Mso for sics- 
plicity w take and to be independent of tij^. WedO so by 
replacing aj^(t)S^(t) by its annual average and by neglecdhg second 
haiisonlc contributions. 

We see that (57) and (59) are Just the equations for the two- 
»de approximation to the mean annual nK>del. Hie only n(» contri- 
butions to these equations are the tc:ms proportional to f)«|(t)Sj|^(t) • 
These terms are called "residuals In this term camss a warming 
of about 2*^C that is not present in the Man annual models The 
residual occurs because the northern hemisphere is more rff lective in 
wintor than in a?OTS/sr (making the seasonal at^lituda lar^) but the 
impact of absorption in sus^r'is greater because the ©mis low in 






th® Bky during winter 0 This wa© not®d by WtMraM nnd M 

Heaebe (1972) ia their ®tuSy ©f tto mmm&l cycln ^th & OOl# 

The seasonal «s^lit^ 2 d© ©f th® sonsl arorage ©wface 
my b« obtained by solving (5d) for Isnd and ocean nassos and .: . 3 

roconstnicting T^(t) ftm (57) ♦ Ths solitude of Tj^(t) for the .. v 

northern hesdspher© is 15«5*C5 the phase lags the solar heating by 
32 days* For reasonable values of the paramters (North and Coakley* 

1978, 1979), !•©•, D taken frota the man annual irsodel; V very nearly 
zero; (^, cossputed heat capacities for a colran of rir end a 

coli^ of ocean missed layer (75 m); « 0*40| the i&odel yields the ; 

correct phase and amplitude for the seasonal ^!^ode» (Note that in ;] 

—2 

this section as in North and Coakley, 1979, we have used Q 340 W m , 

-2 

whereas in Section 2 we used 335 W m .) 

Simple experiments with the solutions show that the asxiplitude 

and phase are rather sensitive to the fraction of land, Tor 

instance, if f, is reduced to 0*20 (southern hemisphere), then the 

L ! 

seasonal solitude Is reduced to about 8®C in rough agreement with j 

the jymetrijsed southern hemisphere data* 

It is possible to parameterize a(x,t) in the model in terms of a 
snow line (or cloud line) attached to an isotherm and there^ produce 
a seT f-consistent feedback scheme* In order to study cliisste change 
or sensitivity, North and Coakley (1979) have considered an ice cap 
whose edge is at the annual "“10^0 isotherm, and a seasionally 

living snen^ line on land whose edge is at the instantaneous^ 0®C • 

isothesna. The resulting energy balance acdel is slightly R>nnHnear \ 


I 
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end the four coefficient# in (51) for T(x,t) are obtained by iterating 
a system similar to (57)-(59). 'fhe resalts Indicate tlmt this model 
has about th® saasa sensitivity to solar constant change# as the 
corresponding mean annual model, the sacll (fe« pcrc^?t) differences 
can be explained in terms of changes in the residuals discussed 


earlier. 
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4. Sensitivity ®ttd the para®et®rl»atien8 

In fcVie precedins sections we introduced sisuple clisaate sodels 
end gave a short survey of solvable easaaplcs. The prissry purpose of 
thees isodals is to guide oar uaderatanding of climt® feedback 
jEechaaisffi®. We saust view with skepticlsa, however, the nuneerical 
results obtained with the 6 l 85 »le isodels. These results are bound to 
be sensitive to the various paraiasteriaations used to link the energy 
fluxes to surface teaperaturcs. In this section we discuss the esoot 
cotissonly used pavasBeterlzations, how they were developed, how reliably 
they reproduce the relationships between the fluxes and the surface 
tes^peratures and, finally, how they influence the fundancntal sensi- 
tivity of the taodel as given by the sensitivity parameter 

To understand how the parameterisatlons affect we note that 
the sensitivity for a global model may be derived directly from the 
condition of global energy balance (Cess, 1976; Llan and Cess, 

1977); 




(64) 


Differentiating (64), we obtain 


dT 

o „ 


I 

o 

(65) 

dQ 

dl 

o 

dH * 

. o 



” ^ dT 

o 



3. 


Obviously 3 is influenced by the sensitivity of the in^-rared radiative 
flux to changes In the surface teciperature, dl^/dT^; it Is also 
Influenced by the sensitivity of the planetary albedo to changes in 
surface temperature, Q dH^/dT^— the albedo-tenperature feedback. 



L®fc# ebT/ioaa# how®v«ir, U iaflMac® of tW mrldloasa Kraassport 

on g . tm mtUicn&l tmmpctt miUcU B ptln^srlly tlns^ngh tu 
o 

inflwaco on tha alfcedo»e®^ei'®tM*‘® fee®>ao5s« «ss.s^l«» th« 

9 urf«c« alfeato is taken to bs a tm&tim of tha local ^usrlaca t®a- 
peratuc®, then edXl be a iwaction of the latitudinal distribution 
of twerature, which in turn la governed by the i»sridies»al trans- 
port. A® a result, the alfeedo-ts^eratur® feedback ia SJebject to 
both the albedo pararaeteri cation end the transport para«teritatlon. 
Kote chat h®**® w« »re considering only the uncartafety In tha 

acneitivity due to a amU ^IgMi ''snE® ^ 

seen, for larger changes a critical solar constant ia r®ni...od beyond 

which only Ice-free or lce-cover«d steady atates eslet. Aa es^haaired 

by t;arren and Schneider (1979) , uncertainties in lead to large 

uncertainties in this critical value. In sisple i^dela seductions 

la 6 tend to be associated tilth lower critical points, iUuatrated 
o 

in Figure!" 13 and 14, but e'/en this correlation could be reversed by 
neglected nonlinear effects, la Section 7 we consider ®.^ll local 
changes in heating and an associated latitude-dependent sensitivity. 
Again, any uncertainty In contributes to a large uncertainty in 
this spatial sensitivity. In the following w^ discuss uBcertsintles 
in the infrared, :o-albedo end ermsport paraiseterizatims, and their 
effect on B^, 


«» tnfrag ffid 

A Bsethod o? psrefiasteriRatiea comonXy «s«d in assergr;;^ 5s«2L«nc« 
cltmtffi tsodeli^^ w®s first «ppli®d by Btidyko (1969) to e! btain th« 
mittod fius in tersa of tha sorface tesperature. E« i)sf;.’?otha8lsed that 
th« way sateoroiogical variables influence the flu* i» «f:!«hibite4 in 
latitudinal end aeaaonal climtological data. To obtain ' the relation- 
ehip he collected meteorological records from a varied cof stations in 
the northern Itecilsphere and with these records he calculested the emitted 
IR flus at the top of the atetasphere* From the reaulta ot'f these cal- 
culations he deduced that the emitted flux is given Isy 

X «> A^ + A^Ag + <Bj 4- , (66) 

where A,, A«, B, ard B, are constants; A^ Is the cloud co«>ver fraction» 
12 12 c 

and T is the su' *ace temperature. Cess (1976) applied (66S) to climato- 
logical records of zonal eean surface tesape*’atures, clouds cover and 
satellite ebeerved emitted fluxes. He found that indeed -<66) provided 
an excellent fit to the cllraatologlcal data. For each 10'® latitude 
cone, with •• 257 W ta”^, A 2 • -91 W vT^, “ 1.63 H sT'^^C and 
Bj “ -0.11 W the paraoeterlzatlon fits the northetrn hemisphere 

data with a maximum ei‘ror of 1®2% of the emitted flux« WIith «» 262 W 
• -81 W • 1.64 W f.^d - -0.09 W the paranr- 

eterization fits the southern haalsphere data with a osaxltaua error of 
2.3S. 

t'Cispite the different climates exhibited by the hemiuspheres, the 
parameterization proved successful at fitting the observatclons with the 
constgjus changing only slightly from one hessisphere to c'he other. Such 
result® argue for the universal ity of the poraseterizatiom. Most 
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r«®artebi® Is that because « 0^ for both hemispheres, Bj « 
ia the only paraseter, according to J65), that would ef fact 0^^ and it Is 
the eaif® 1.6 W a”"*C for both Healspheres. The s<mwhat larger value 
of B " 2.09 W used in earlier sections for illustrative purposes 

srisea fro® applying I • A 4- BY to the sonal clieatological data and 
ignoring cloud cover. Also ignoring cloud cover, Oerlcaans end Van 
dsn Dool (1978) obtained 2.23 W 

Sellers (1969) developed anc.her parataetarltation for the eottted 
flux. It la nonlinear in sorf*>'v teaperafure North (1975b), however, 
showed that the paraoeterlsatlon «iicn linearized reseasbles the 3udyko 
paraeseterleation with constants that oiXfer only slightly frora those 
given by Budyko. 

That the emitted IR flux is so simply related to surface temperatures 
seesis miraculouo when we consider that the earth’s surface and each seg*" 
(sent of the atmosphere contributes to the emitted flux. The contribution 


made by a Bcs®ant of the attaosphere depends on Ito temperature and con- 
centration of emitters. Water vapor and clouds are the saajor emitters. 

As clouds and water vapor are confined t *, the troposphere and ft« tropo- 
spheric temperatures are sufficiently high, most of the I.R flux emanating 
from the atmosphere la emitted by the troposphere. Thu* My feedbacks 
that link surface temperature to tropospheric temperature profiles, 
humidity or cloud cover are likely to influence th- sensitivity of the 
outgoing flu;'? to changes in surface temperature. 

To llivMstrate how various feedbacks might affect 0^ we take dl^/dT^ 


to be given by (Cosklay, .197?) 


dl„ 

dX. 


31, 


dr 


"o "o *'■ "o 


(61) 
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In (67) <SA /(!T represents the rate of cloud esEtount change with surface 

CO 

temperature chmge,> tte clou4 raorat feedback; df MX tha cloud top 
tcoperatur® feedback; dF/dT^ the tropospheric lapse rat® feedback. The 
partial derivatives in (67) are evaluated uihlle keeping the riasalning 
variables fl^ed. 

Feedbacks with the largest potential influence on dZ^/dT^ are 
those related to cloud coveu'. Cess and Iteanathen (1978) report that 

m»2 

-91 M o < 31^/9A^ < -34 W a , depending on how the vertical profile 

of cloud cover changes as the total aisount of cloud cover changes. We 

•*1 

note then that a cloud as»unt feedback dA^/dT " 0.02*C , would con- 

G O 

tribute as asuch as 1.8 W to dl/dT . Tliat lo. Its contribution 

0 o 

vould be K8 large m the v.'lue deduced for dl^/dT £ro» climatological 
^ o o 

records* Clearly* cloud amount feedbacks could greatly influence the 
sensitivity of the climate* 

Another cloud related feedback is the cloud top temperature feed- 
back. Calculations made with global averaged vertical column models of 

the earth’s atmosphere (Ramanathan and Co^fkley* 1978) indicate that for 

-2 —1 

fixed cloud top temperatures* 31 /3T « 2.16 - 1*75 A or 1#29 W m 

o O C " 

for “ 0.r>a while for fixed cloud top altitude and constant tropospheric 
lapae rate - 1.0) (dT^/dX^) - 2.16 0.19 A^ 

or 2.26 M for A - 0.5. Thus, dl /dT could change by 1.0 W 

C DO 

if the cloud tops Utsnge from Mintaining constant cloud top temperatures 
to maintainiuz constant cloud top altitudes* Hw the clouds might change 
during a change* however* remains a mystery* 

The prec&dii\g values for dl^/dl^ were corsiputed assuming that the 
maintaln^i' eon^tmt ptoBlae of relative, lnjsuidity. \: 
the. tropospheric temperature "Increases tka concentration of water vapor ' 
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'I 

■ ■■ / -'i 

Is 88su®sd to incteffisa. Th« a8»uss|>tlon Is baaed on the «hs«rvatie» that | 

the aatth’a stesosphete appears to conserve relative Iass8i#Ity (Kaaabe and 
Metherald, 1§67), Global circulation w>del« of the earl's atssoaphere 
also eeea to conserv® relative hisaldlty (Maaabe and Wetlarald» 1975s 

Hatharald and Hansba, 1975). If, on the other hand, tha ;«raosphere j 

••2 “1 ' 

©aintsined profiles of constant absolute htauidity dl^/dXjj ^ 3.7 W B *C i 

under the conditions of fised cloud top altitude (nassnathm and Coakley. | 

1978). We note then that without the observed ooisture feedback affecting 
the opacity of the atsaosphere, the emitted radlatisu esshil^fesa the sensl- ; j 

tlvlty of a blackbody at the earth's equivalent temperature • 254 K. 

Finally, as with the ©oisture and the cloud cover feedbacks, lapse 
rate feedbacks ar<* also expected to Influence dl^/dT^. Raaimathan (1977) 
noted that Ip the 0-2% solar constant change experiments peafbrmed with a 
GCM (Wetherald and Kanabe, 1975) dl/dT with fixed cloud top altitude , 

ranged from 2.4 1.45 near the equator to 1.7 ~ 0.81 near the pole. 

He attributed the range to differences In lapse rate change&i Near the 
eqvuitor moist adiabatic adjustment governed the lapse rate ohanges while 
near the pole a mixture of radiative and advectlve processeas governed 

the lapse rata changes. | 

Obviously, the feedbacks mentioned thus far could have® profound , 

influence on dl^/dX^ and thereby on 6^. How these f eedtacke work and 

what other feedbacks might affect the climate sensitivity rspresents the * ; 

focus of much ongoing research. Given the range of possibf-S.tlea ■ r 

tUustrated thus far, it Is difficult to specify a probable range for | 

dl /dt . Barring cloud amount changes, hcswcver, we might aeapt Che | 

o o ■ ''. 1 ? 

awl 

range from 1,3 W eT ^C obtained with fixed relative huaidiy, lapse ; ' 

-r»2 ■ I 

nste and cloud top temperature to 3.1 W m *C obtained wife fixed ’'I 

• ' ■ ':s 

V ; ' ■ '■ ■ ^ ^ 


relative husaidity^ cloud top altitude and laoiet adiabatic adjustsicnt 
OKhibitcd in th® tropical region by the GCM. 

In accepting this range* ^ should note that contrary- to the finding 
by budyko and Cess of a universal constant for erperis^nts vith a 
CCS! Indicate tliat dl/dT changes with latitude. Oddly enough, the Budyko 
parasaateritation fits the tonal clieatology of the GC21 about as wall as 
it does the earth’s (Coakley and Wielicki, 1979). The apparent discrepancy 
is resolved when we recognise that the feedbacks that influence the clioate 
change exhibited by the CCH are not revealed by applying (66) to the tonal 
fields. In a sltailar way, we oight expect that application of (66) to 
the earth’s fields would also miss the feedbacks that could significantly 
affect climte change, 

b. Albedo parameterization 

As the earth cools we expect the extent of penaanent ice and snow 
cover to increase, thereby Increasing the earth's albedo. To simulate this 
effect, Budyko (1969) proposed the siople mechanism Introduced in Section 2. 
Poleward of the -10“C mean annual isotherm the surface is to be covered 
with ice; equatorward, it la to be ice-free. The -10*C Isotherm represents 
the boujidary of permanent ice and snow cover for the northern hemisphere, 
Budyko assigned albedos of 0.62 for ice-covered regions and 0.32 for ice- 
free regions. Hence the albedo change is 0.3 when the surface changes 
from bting ice-covered to ice-free. As la indicated in Figures 13 and 14, 
such a large change in the albedo at the ice line makts the model highly 
sensitive to aolor constant changes and causes it to produce the completely 
ice-covered solutioti when the solar constant is only slightly reduced. 

him and Cess (1977) noted, however, that because of the zenith 
angle*n',ependent reflectivities of clouds and surfaces, albedos of the 
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ice-free regions will be higli at high latitudes due to large tenlth angles 
and low at low laticudea. As a result, th® chisnga in albedo as the 
surface changes froa ice-free to ice-covered conditions will be large at 
low latieudea but saall at high latitudes where the ice cover changes 
taka place. Allowing for the senith engle-dspendent reflectivities referred 
to by Lian and Cess, Coakley (1979) finds that the change in albedo at 
the ice line is reduced to 0.15, for which QdH^/dT^ - 0.4 W under 

the current climatic conditions. 

Instead of an ice-line. Sellers (1969) allowed for the change in 
albedo by taking 

+ c'T T < 283.16 

( 68 ) 


a 


{ b + 


2.55 


T i 283.16 


From albedos and temperatures in the same latitude zones of the two hemi- 
spheres he deduced that c - -0.009 He then made b a function of 

latitude so that (68) matched i^bserved albedos. 

Because of the zenith angle-dependent reflectivities noted by Lian 
and Cess, however, we would expect c to be a function of latitur’e and not 
a constant. From climatological observations of zonal albedos and surface 
temperatures, they deduce that c ranges from 0.0 for latitudes equator- 
ward of 40*N [near the position of the -»-10*C Isotherm and In agreement 
(68)) to -0.0U5 K-' « 85-H. U.lng the l.tUudlo.I-depenient c. 
thc> obtain QdH^/dT^ " 0.3 H m *c ^ while using the eawe model but with 
c » -0.009 K ^ at all latitudes, they obtain QdH^/dT^ « 1.0 W 
This difference taore than doubles the aensltlvlty of the model. 

As with the emitted radiation, water vapor strongly Influences the 
absorbed solar radiative flux. For fixed relative humidity, absorption 
by vufater vapor increases as the tropospheric temperature increases. 





Coakley end Wielickl (1979) ©stimte that for fixed relative humidity the 

change in absorption by wtor vapor contributes 0,2 ^ to QdH /dX • 

o o 

Of the solar radlatioa reflected by the earth, clouds reflect between 
70 and 80%, Wa would expect, therefore, that cloud avmmt feedbacks would 
strongly affect the sensltivley of the absorbed solar radiati^'e flux to 
changes in surface temperature. For global average conditions (Cess, 1976) 

Q3H^/8A^ » QCoig - a^) , (69) 

I'lhere a » 0.18 is the albedo for cloud-free regions and a « 0.43 is the 
0 c 

••2 

albedo for rloud-covered regions# Henc©» Q3H /3A « -^5 W © and thus 

o c 

*•1 

dA^/dX ■ 0,01*C is sufficient to swamp the feedbacks examined so far. 

Clearly, the cloud amount feedback could strongly influence the sen- 
sitivity of the albedo to changes in surface temperature just as it could 
strongly influence the sensitivity of the emitted IR flux. We note, however, 
that the change in the absorbed solar radiative flux caused by 'a change in 
cloud cover is somewhat compensated for by the change in the emitted IR 
flux. That is the term In the denominator of (65) contributed by a cloud 

amount feedback, (3l /8A - Q3H /3A )dA /dX is smaller than either 
o c o c c o 

31 /3A (dA /dX ) or Q3H /3A (dA /dl ). It, in fact, may be negligible 

(Cesa, 1976)# Its tnagnitude, however, is the focus of considerable debate 

(Cess and Raraanathau, 1978; Ohring and Clapp, 1980; Hartmann and Short, 1980)# 

Barring again cloud amount feedbacks, we might accept a range from 
-2 -1 

QdH^/dT^ w 0.2 W m *"C for fixed relative humidity and no ice-*albedo 
-2 -1 

feedback to 0.6 W m ®C for fixed relative huHiidity with ice -albedo 
feedback. The larger values obtained with the Budyko end Sellers parato- 
eterizations appear unwarranted in view of the zenith angle^^dependent 
reflectivitiy cox’tectious which they neglect# 
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c. Trans^ , ■ 

In the two-«ode approximation the sensitivity of sliaple taodels was 

shown in Section 2 to be unaffected by the model used for the meridional 
transport. Held and Suarea (1974). Horth (1975b) and Un (1978) have also 
shown that various diffusion parameterisations that are nonlinear in the 
meridional temperature gradient, such as that suggested by Stone (1973). 
also have relatively little effect on sensitivity. If a third mode is 
added, then corrections to of order 1C% are to be eKpected. 

Allowing for a Hadley cell along the lines suggest ‘d by Undzen ano Farrell 
(1977), for example, might be expected to make such a correction. 

Remember that in the simple models the transport only affects 6^ 
indirectly. TtiC transport enters the energy balance as the divergence 
of a heat flux and integrals over the globe to obtain the global energy 
budget cause the transport tern [-R^DV^t or y(T-V 1 to disappear from 
the global energy balance (64), and consequently because dl^/dT^ is 
usually taken to be constant, the transport can Influence the sensitivity 
only through its impact on the albedo temperature feedback, QdH^/dT^. 
Without albedo temperature feedback, the meridional transport could have 
no Influence on the sensitivity of these simple models. 

A realistic model for the transport has yet to be constructed. In 
an attempt to check the diffusion vaechanism, Lorenz (1979) has recently 

a study of sensible heat transport in the real atmosphere. He finds 
that diffusion parameterlzations may work for the largest scales (Legendre 
index n - 0,1,2) but surely fall for smaller scales (flux and temperature 

gradient are improperly correlated). 

Despite the osmtprcsence of diffusive transport models, the results 

of the GOT referred to earlier suggest that Its transport remains constant 
independent of surface temperature and regardless of external forcing 
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*“ to the l„s» clwaj,, ^ r th. c® i., 

^ the '«' ‘01„, 

“ t. the .,„tUh«„.°eZr fT 

p«v.te c«.„„tee«e„> -0»1 <H«..be. 

eft... Rooeht etudiee „. a„ oi“«lhrl„ 

»Mm to bear thla cot. ™e<Hae,.,ij,j arolel by Held (uyg) 

Wth the rangea obtained for B . dl /dy , , 

«nd for QdH /{Jt « , o o» < B < 3 .x w 

o^o> 0,2 < QdH /dr < 0.6 w »-2.,o ’ 

“•«'c < B. < 3.4 "c ii.h a ° ' "“ “ '“"SO 

o c. although difficult to .oMf, ce, 

“^t» the actual oeualtlnty of the cllaate f' 

•ooottlvltlea of acc clloate t»del.-e„ccsy b^ 

thoee „lth the l.rg. albedo-teaperature feedbalar' 

“Oola (Ramenuthau and Coaldey, 157s) ’ "'“‘“'"-oo-vectlve 

•oOele (Lieu and Ceaa. 1977, “ “‘ttolatlon cllnat, 

'ooPonae of the cllnate to 

'oe^t Of unicanic eruptl„„ «ch„elder «,d “ ‘ 

Ranaen „ «■"»• t,75,-Ha.a and Schneider, 

the effect of ,t leapt cme ' 'o«e«hered. houever, that 

Jtloned before, uhat cloud cover change, ufn 17 ^ ^ 

thange. and ho„ they «in obloate 

research. ® subject of current 
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5. Sensitivity to changes in orbital parasisters 

Over periods of 10^-10^ years, the Earth’s orbit about the sun 
changes, Miloakoviteh (1941) aaong others argtied that these changes 
cause variations in the asfflunt sad distribution of solar radiation 
received by the Barth, thereby influencing the clitaate. Milankovitch 
suggested that the orbital changes force the advance and retreat of 
glaciers. Indeed, evidence that the climte is so forced has been 
found in geological records. By examining sediments taken from deep 
ocean cores. Hays ct al, (1976) have shown that the principal periods 
of climatic variation (100,000, 42,000 and 23,000 years) match those 
for changes in the eccentricity (105,000 years), the obliquity 
(41,000 years) and the longitude of perihelion (23,000 and 
19,000 years). Many have used simple energy balance models in 

efforts to explain the glacial cycles as the result of orbital 

changes e 

How energy balance models respond to orbital perturbations is 
understood through the effect of the perturbations on the forclng- 
the incident solar radiation. To start, we note that 'fractional 
changes in global average m an annual incident aolar radiation are 
approxiMtely given by 1/2 Ae^ (Berger, 1978), where e is the eccen- 
tricity of the earth's orbit. Because the eccentricity has always 
bftv.i small (e 0.07), the resulting changes in the solar constant 
have always been less than Q.2%. A 0.2S change in solar constant 
causes, according to the energy balance taodels, about a 0.4“C change 
in the global wan surface temperature. Such a change is an order 





of smaller those that mnm to have taken place (Hays 

et 197S)« Ctelng to th©ir insignificant magnitude^ changes in 

ener^ balanc® models brought about by variations in the eccentricity 
will be mglected in the follwlng discussion* 

Ifore important than changes in the m^&n annual solar radiation, 
however, ®re the variations in its latitudinal &nd seasonal distri- 
bution# thea® wriation© result from changes in the ©bli<?ui*wy and 
longitude of perihelion* The seasonal variations are clearly repre- 
sented by changes in the incident solar radiation for summer and 
winter caloric half-years, as defined by Milankovltch (1941) (Berger, 
197£). may derive these changes from the expression for S(x,t), 

S(x,t) - S^(t> + S^(t) F^(x) + S^(t) (70) 

where, to first ordar in e, 

S^(t) * 1 + 2e cos (2iTt - II) , 

Sj^(t) "» S^Ccos 2TTt 2e sin II sin 2Trt) 

and ' 

S^it) - ^ cos(27Tt - 11)1 . 

Here, t ** 0 at the northern hei&aisphe re winter solstice; II Is the 
longi :ud® of perihelion is^asured from the longitude of the winter 
solstice, and In the derivation of S^Ct), we have neglected the first 
hsr^>Bic terms proportional to o that arise from the second harmonic, 
which res^islts fro^ the nun crossing the equator twice each year* 

Xlie deriv^vtion of the tetM In (70) Is outlined by North and Coakley 
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(1979)* Averaging S(K|t) over the suisiser and tdlnter half years^ 
we obtain 

and 

®w " ^ S^P2(x) + ^ ainn [1 + , 

Changea in the Incident radiation are thus given by 

ASg - - f as^P^(x) + AS 2 P 2 <x) - ^ A(e sin Jl) [1 + S^PjCx)] (72) 

and 

“ ? ^S^PjCx) •*• AS 2 P 2 (x) + ^ A(e sin H) [1 + S 2 P 2 (x)] . 

The above expressions are the three isode representations of the 
Mllankovitch forcing as it has traditionally been calculated (Mliankovitch, 
1941; Vemekar, 1971; Berger^ 1978). 

From (72) we see that the change in the mean annual distribution 
of incident radiation is given by 

AS(x) « I (ASg + AS^) « AS 2 P 2 (x) . (73) 

S 2 la easily shown to depend only on the. obliquity, 6^. Because the 

range of obliquities remains small, 22.10® “ <5^ ~ 24.51® (Vemckar, 

1971), we may take Ac m 0.0115 A6 , where Ai*' is in degrees of arc. 

^ f o 

Currently^ 6 « 23.45*^, but 25,000 years ago 6 was about 22.2®. 

o o 

This change in obliquity caused there to be less radiatlan Incident 
at the poles and more Incident at the equator. The change Is thought 
to have spurred the glacial Bsaxltmira 18,000 years ago. Inserting the 
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altered distribution into a mean annual Siodel, however, gives only a 
1—2* equatorward shift in the latitude of the ice line and a wanaing 
in low latitudes (Budyko, 1969; Sellers, 1970; Saltzman and Vemekar, 
1971; Coakley, 3979). These results clearly disagree with the ice 
line shift of approximately 15® end a cooling of several degrees at 
all latitudes deduced for the glacial maximum (CLIMAP, 1976). 

The small changes in the mean annual distribution of incident 
radiation and the disappointingly small response of mean annual 
models has led some to suggest that seascnal variations need to be 
considered when estimating the effect of Milankovltch forcing. 
Obviously, (72) indicates that variations for each season will be 
larger than those for the annual mean. Results of seasonal models 
that allow for these variatljns, however, fail to differ signifi- 
cantly from those of the mean annual models (Suarez and Held, 1976, 
1979; North and Coakley, 19/9). to first order in e, is easily 

shown to depend only on 6^, and thus, like ^ 2 * find “ —0.016 A6 
Because S^ contributes only a small lncreiL‘iit to and X 2 » ASj^ 
raises the response of a seasonal model only slightly over that of a 
mean annual modal. Likewise, by considering as in Section 3 the laean 
annual residuals contributed by the terms propoi -lonal to e in (72), 
we find that for the last 25,000 years these terns contribute at ciosc 

variations of the order of 0.1*C in T and T-. Also, if we allow 

o / 

for differences between the annual mean temperature of the two hemi- 
spheres, the terras that couple the seasonal variation of albedo with 
changes in the longitude of perihelion contribute at most 0.2“C to 
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thi?' mm ccstgponent of ThoSt aroual not 

6®i.,*:DnaX i!^d€^ls giv@ tmpomm that wuld support the Hllankovitch 

theory# 

Of courses p as is^ntioaed in ihe two previous se^ctions* sls^licity 
and questionable para^teriisations my be the reason that energy 
balance /ail to give large responses. Perhaps if more phys** 

leal procesfii^a were ?^deled, more feedbacks Includedt laore realism 
added» the models would become naore sensitive. 

Sow additional feature*# have been recently added. Pollard 
(1978) added Weertmar.’s glacier mdel (Weerman* 1975) to a ^liapie 
three-'xaod'R seasonal saodel. Becatase it allows for the ablation of 
snow and ice by Incident solar radiation, the glacier model contributes 
far t?ote to the sensitivity of the seasonal rsodel than to that of the 
mean annual isodel. The enhanced sensitivity is Indicated by the 
iscnsltivity of the Bscdel tc changes in perihelion rhown in 
Figure 15. Keverthelese, if we scale the results gshowtt in Figure 15, 
we find that the TOdel, even with enhanced sensitivity, would predict 
only a 4-5^ shift In the ice line for the obliquity change connected 
with the ia,CD0 TBP glacial maximum. The addition of a glacier model 
by itself hi>H so far proven insufficient to obtain the desired 
results. 

Another factor to consider is the zonal asyts^etry of ice ag® 
climateis. Hotu of The Ridels mentioned so far have been «;onslly 
syiraetric. Reconf--tructions of t\m glacial indicate, 

however, that in tha mtthem homisphtr^ ice extended mch further .. 
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««iu®torMaKd ovar than ©var oc«ana (mmP. 1976). mttmti 

and Short (1979) hav® «ti«ad that th« glacial ®aai»a and its evident 
mymmtty my hav® bsen ths to a aaywstrlc wavs in 

ths auriace Using sn anatg^ baianca that eontalna 

a prsacribseS mymmtry in ths surface te«|»®rsture» they obtain an 
lce-ag«»llh« response when they Insert tMSonable «apSiwdes for the 
asywtrlc t«m. With this additional cowownt. however, the «odsl 
appears to be no »r« sensitive to changes in the lncl*^nt radiation. 

In Figure 16 the global seen tes^erature la plotted as a function of 
solar constant for different «« the tonally «y«®etrlc 

coaponent. Th® nsacly Identical slopea of the curves isdlcate a 
co»on sensitivity. From this result we «*P«ct that allcwlns for 
the tt6y®a®try would not greatly enhance the sensitivity ot the model 
to orbital changes. In f.ct. by allowing temperatures ever continents 
to differ from those over vcenns. Sunret and Held (1976# 1979) 
incorporated a degree of a.y^etry In their model, nve results of 
their model, however, are not significantly different from those of 
ronally symmetric models. 

Still, failures to euppcrt the Milankovltch theory may only 
reflect the inadequacies of the models. Moat models ttet have been 
used to test the theory have Included only one feedback that could 
amplify the response-^the olbedO“te»perature fecd^ • 
with these model®, however, leads ua to search for eddStlonaX feed- 
backs that either by themselves or working together cnrld aapUfy the 
response (Coos and yronka. 1979). 


6* Stfibility theory 

Wm A m thoM dincmm4 in Section 2^ h&m hmn 

aoived md tseveroi @ointion» ere for the ©am eKterael condi- 

tioner m ©hould stady the etebllity af each aolwtion. Schneider 
end 6eI-Chea (1973) veto the firat to @tndy the effect of vsrioue pertur- 
bations on a Sellets-type md&l at fixed eolar constant* They found 
that the present climate va© ©table? but that if a large enov ,h cold 
perturbation (^20®C) was applied at all latitudda, the solution 
failed to return to the present climate but Instead pltmged to the 
ice-covered earth solution. Subsequently^ linear etabillty results 
were obtained for the Budyko KK?del by Held and Suarer (1974) and Su 
and Haieh (1976); other analyses for diffusive saodels were performed 
by Horth (1975a), Ghil (1976), Fredericksen (1976), and Drasin and 
Grlffei (1977). The ^st general linear stability analysis for a 
wide class of models was perforrod by Cahalan and North (1979). An 
approach for a finite sm;>litude stability analysis was suggested by 
Ghil (19/6) and a coiSEplcte study along the same lines based upon a 
variational principle was given by North (1979). Incorporating 

ideas taken fi'om all of the above-mentioned studies, we shall now 
discuses the stability of the model solutions. 

In performing a stability analysis, we study the tte^-dependent 
respon^ie of the to perturbations from its eqaiiibrlura state. 

We shell restrict the discuission here to t&ean arimaal laoJels. A cli- 
tsatic state not in equilibrium wjst have a heat storage term added 
to TOergy balance equation, C ST(x,t)/3t, where C is the heat 
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capaeity par uatt area for a eoltm of the ®fflrth-at®o»phare eyatea, 

As santlosaiiii la Sactioa the vast difference in C over lead and 
oeesn leads to Bcm error associated vith ei^ly taking a sonally 
averaged C» but m shall ipiore thie co®g>licatlon la the present 
section. It is easily seen that C nerely ecalea the tiosa in our 
present discussion and since ve are only asking whether or not a 
given solution is stable ■ regardless of tfjw scale. *^e taay set C 
equal to unity. 

•• Linear stability of global models 

To illustrate the concept of etebllity. consider a oodel with no 
ice cap feedback. Its tlae behavior may be otwdied by adding the 
storage term to (8): 

3T (t) 

+ A -f BT^(t) - Q(1 - Op) - QH^ . (74) 

It la easily shown that any deviation from the equilibrium solution 
given by (8) decays back to equlllbrlitm exponentially .with a time 
constant equal to 1/B (times the heat capacity per unit area C, which 
we have set to unity). The outgoing radiation therefore damps the 
solution back to equilibrium. For a column of attaosphere only, this 
relaxation time Is about 58 days. Larger values of B (Increased 
negative feedback) lead to Increased damping. With no albedo feed- 
back the laodel is linear, and it has only one solution that is always 
stable. 
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The nest b»s!s 1 of Interest slXowa eoae tflwsperature dependence in 
the planetary albedo; It i» the sero-diaeasional Ice cap »odel dis- 
cussed earlier (cf»» B<j. 13) • In this case 

<t) 

« QHolXg<V^ • 

and th® ssodsl tay be linearised for aaall deviations from a ateedy 
state Bolutioa, say T^°. Then If ve let 

T (t) - T ° + 6T(t) , (76) 

o o 

we may write 

dH 

6T(t) * B 6T<t) - Q «T(t) , (77) 

where we have made use of the steady state condition given by 


A-hBT^ - Q»o<0 


(78) 


Taking a derivative along the curve given by (78), we obtain 

!< . „ < 

dQ "» ldT° “0 


which In turn my be used to simplify (77) giving 




6T(t) ■ 0 , 


(79) 


(80) 


Equation (80) establishes a theorec which may be generalized to a very 
large class of clissate models; If the steady state solution is on a 
branch with dX^/dQ positive, the solution will be stable; If this la 
not ao, the positive albedo-tetaperature feedback, Is larger 
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them the XE radiative dsa^ing B end the eolution will h# unstable. 
To sm^rise for the si^llfied global sadel: 


»£ > 0 stability 


.«. < Q ixkstability 


Xt is also fairly easy to exasaine qualitatively the ^msponse of 
the BlBsplifled global eodel to large (nonlinear) perturbad-ons. We 
present the analysis here as it will prove useful In the oore cosq>ll- 
cated exas^les to follow. Our itsethod will be to construct: a potential 
function (sometlE^s called the Lyapunov function, cf*. WorJh et al. . 
1979) that will completely describe the laodel and Its behavior away 
from steady state. 

Consider the function 


F(T„) - AT„ + i - <M„(T^) . 


where 


y ^ i: 


H (T^') dT^’ , 
o o o 


is detertsined from the tempera ture-dependent co~albedo shorn in 
Figure 1. A graph of the potential F(T^) is shown !.« Kigu® 17. 

Kote that the axeretsa of P (points where dF/dt^ • 0) are gh/en by the 
energy balance equation (78) ; therefore, these extrema cor%spond to 
the roots I, I'i, and III of Figure 1. Note also that in tanas of 
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P<T ) (75) becosaes 

o 

T - -dF/dT . 

o o 

Furtharmora, • . 

dt d?^ 'o o' ; 

Thus, the tide derivative of is proportional to the slope of F at 
T and F always decreases with tine. The Implication is that the 

shallow local minimum 1 (the present climate) Is stable, the maximum II | 

is unstable, and the ice-covered planet III is stable. Figure 17 j 

makes it entirely plausible that a uniform -20*C perturbation applied j 

to the latitude-dependent models would probably take the solution | 

"over the hill" and into the ice-covered planet valley, as was die- | 

covered numerically by Schneider and Gal-Chen (1973). | 

1 

Note that the "slope-stability" theorem (81) and the existence | 

and qualitative form of a potential function (82) are independent | 

of the exact parameterlzations used for the albedo and IR radiative j 

flui:. As noted in Sections 2, 3 and 4, much remains to be done in | 

the area of parataaterleatlon theory so that results independent of j 

details or numerical values of constants are of special interest. . | 

In the remainder of this section we will Illustrate how the slope | 

stability theorem and the potential function may be generalized to V| 

'I 

the class of one-dimensional models. We begin with the linear ^ 

I 

stability analysis. 1 
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c» Linaar stsbilitv of one-t^issnsAonal gsodals 

The ti8®-dependent energy balance equation la pl'^n by 

0(s)(l - x^) l«,t) -b l<x,t) » Q S(x) a(x,Xg) 

( 86 ) 

and the ice line condition by 

■ 1 . (87) 

9 8 


In theoe equations ve have used I • A + BT rather than T as the 
dependent variable* The tlwe In 3l/3t has been scaled by C/B, whert 
C Is the heat capacity per unit area. Expanding X in terms of the 
eigenfunctions of the diffusion operator, ®^y wse the ortho- 

gonality condition (47) to determine the coefficients. From (86) 
ve obtain 


1 + HI - Q h (X ) , 

n n n n 8 


( 88 ) 


where we have defined expau^ion coefficients for 1, 


t'n< 


In(t) " j f_(x) I (x,t) dx , 


(89) 


anci for the solar heating, 


n s 


f f 


h (x ) “ f (s) S(x) a(x,x ) dx , 


(90) 


and 1 ■ n +1-1. 

n n 

To linearize about the present climate let 


I (t) « I ^ ^ 6l it) , 

n ' n n ’ 


( 91 ) 
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m& 

« (t) “ r + 6 js (t) • (92) 

B 0 0 

with the state given by 


a I ® • Q h (:c ) , 

n n ^ o' * 


the ice line condition by 


I! C • 


(93) 


(94) 


These latter cquatiOKS define the cquillbriua state denoted by ) 

vhose stability is to be tested. Linearizing (88) we obtain 


6l + £ 61 “ Q b’ (x ) 5x . 

n n n ’no o 


(95) 


The perturbed ice line 6x can be related to the 61^^ through (94): 

) (5 • 

In (96) and in what follows we suppress the argument x^ in f^, f^, 

h and h*. Substituting (96) in (95) leads to an infinite set of 
n n 

simultaneous » hosac^geneous, linear, first-order .differential equations 

for the 6I (t). By substituting 6I (t) - 6I e"^^, we way find the 
n n n 

set of values of ^ for which the equations are satisfied. This turns 
cut to be tb«> eigenvalue problem 


\ M 6I •» X6l , 
Zrf nit m n 


( 97 ) 
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4^*. V^ >• •. ji. •' ? *■ .*> ■'"*•■ 


. (99a) 

slbedo to be a ^tepfunction In x* For 


. . VhiB' 


;am- \;h^ <>u) can be differentiated giving a delta function 


; -tbs' vOitegr^'J/ .uign * lead'i-^a; the resxilt 

h^' - La S (x^) » C99b 

vhere > 0 is tne discofstiiui_..y In a(x,x^) at x " x^. 

■ I'he Btablli^y vf • the ays ten with steady state i® 

detfemincd by the si^a of the eigenvalues X. Because is real 
and Bycaietrlc, all eigenvalues are real and bounded from below. If 
the lowest eigenvalue is negative, the system Is unstable; 6l grows 
exponentially with time. If the lowest eigenvalue is positive, the 
solution is stable; 6X decays exponentially with time. By casting 
the eigenvalue problem into a different form, we can determine the 
' Ig'^. of t loot* 

To determine the sign, we rearrange (97. ...id \ise (98) to obtain 


(£, - X) 61 - - Yf / f • 

' ^ n n ffi Bi 


( 100 ) 


Dividing this expression by ~ multiplying by and suasaing 


over civ. "vin 


fii 


n V ■; 
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This relation is a trancendcntal equation that Is satisfied for certain 
discrete values of X, the stability eigenvalues. By further rearra'’qe- 
Bjent ve arrive at the sign of the lowest eigenvalue. With 


(99a), (101) becoKOS 

2 


l(r4 

n \ n 


I ® f * 

n rt 

QbS(K^) 


which with (93) becoaies 




) 


h f ' 
n n 


) 


( 102 ) 


Since 


dx 


dx 


[> ? ¥'1 ■ 


h f ' 
n n 

i 

n 


« 1 ^ T - I 

Q dx n f 

^ o n 


(103) 


Substituting (103) and (99b) into (102), we obtain 

bS(x ) f ^ 


dx 

o 


I 


s n 


£ (£ - X) 
n n 


(104) 


As a function of X the right-hand side of (104) has a zero at X " 0 

and a sequence of poles with positive residue at the points ^ 

(£ > l)a This is sufficient for us to nake the schematic plot shown 

n 

in Figure 18. If the horisontal line corresponding to the constant 


dQ/dx^ is also plotted on the graph, the roots X 


(j) 


are at the 


Intersections of these curves. Clearly If dQ/dx is positive, all 

0 

roots are positive and the solution Is stable. If dQ/dx^ Is negative, 
the lowest root becomes negative, and the solution becomes unstable* 




...0 . V . . , 



These results represent th® slope—stablllty theorem fee the siiaple 
cme-dlSKnslonal g»del8. 

The proof sketched above can be generalised and wda rigorous 
for ail sodels that have positive Green's functions (c&» Eq. 49) 
(Cshalan and North, 1979). Note that the slope-stability theoreia 
holds only for north-south symetric solutions and for sKJdels with 
lce~caps that follow isotherm, i.e., the co~albedo tsaybe written 
as a(x,x^). Models possessing an albedo with the funcl^ianal form 
a(T) may have solutions that are unsyraastrical (Drazin md Grlffel, 
1977); in such cases the proof of the theorem falls. Tie proof of 
the theorem requires the existence of the single Index x that picks 
out which of the several solutions one is examining for a given Q. 

It should be noted that the slope-stability theorem for the 

one-dimensional models differs slightly frem that for the global model 

(81). The sign of dT^/dQ and dx^/dQ may differ near a cisp at 

which dx^/dQ « 0. To see how tjiis may happen, we dlffer«tiate 

(93) for n « 0,(A + B T ^ Qh ) to obtain 

o o 

BdT^ dh^ dx^ 

dQ “ Q • (105) 

o 

Both h^ and dh^/dx^ arc positive, Therefore, dT^/dQ readhs positive 
near a cusp where dx^/dQ is small and negative. In that CiBse the 
solution is unstable even though u pc'itive. On the :<ther hand, 
near a bifurcation the magnitudes of dx^/dQ end dT /dQ are large and 
of the same sign. So the positivity of is always necesary for 
stability, and if we exclude cusps then it is sufficient » well. 
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<5* Potential functional for one~dlt3ensiort»l g»del8 

We consider now the finite ajBplitude stability analysis of the 
one-diraensional a(T) saodels. Our aim Is to construct a potential 
function analogoas to (83). In the rero-dliaen-iional case the 
potential (Lyapunov function)* F(T^), was a function of only one 
variable; its local extrema gave the steady state solutions to the 
energy balance equation. 1^. the solution was perturbed, changed 
in time according to the negative slope of F(T^), which Itself 
continually decreased in time until a (steady state) extremum was 
reached. Clearly a local minimum corresponded to a stable solution 
while a local OKOtlrauro corresponded to an unstable solution. 

The aiialogous procedure for one-dimensional models is more 
complicated. The potential must be a function of T(x) at each local 
point X. That is, it must be a functional, F[T]. The nature of this 
kind of mathematical object is most easily described by example. We 
shall, therefore, present a functional and show that it has the desired 
properties (North et ed., 1979); a similar method was auggested by 
Ghll (1976). We take the functional to be given by 

F[T] - / dx [| D(1 - x^) T^^ + R(T) - QS(x) C(T)] , (106) 

where T^ “ dT/dx , 

R(T) “ I(T') dT’ , (107) 

and 

C(T) “ a(T') dT’ . (108) 
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It le relatively straight forward Qorth e^ al » ^ 1979) to that If 

6F is to vanish (local esEtremm) for an arbitrary varlattem 6T(x) in 
the temperature field* the corresponding energy balance eoiuatlon must 
be satisfied (steady state solution)* North et al. (1979) /also show 
that these solutions are stable if the e^t iBum is a lo '.al E 5 sd.ainiuffl 
and cnf^.table if the ejctremum is either a saddle point or a :iocal 
maximum* 

Here we shall illustrate the procedure for the modo^ Alust used 
to prove the slope— stability cheorem* ITie co— albedo* howeve,*r* will 
be given by 

a(I) - e(I - 1^) + a^e(I^ - 1) . jl 09 ) 

where 6(Z) is the unit step function* 0 •« 0 for Z < 0* and «» 1 for 
Z — 0. Note that the step function albedo (109) is the one^cas^. where 
models having the coalbedo functional dependence a(T> and mozsdels hashing 
the functional dependence a(x*x^) are Identical. Xt is easiver to visualize 

O 

the functional F(A + BTJ - F[I) In spectral form, I » Z • Xr 

tpectval form F may he thought of as a function of the variaibJes I ,1-,I.,.. 

o 2 A 

where for s^nsxmetric models we use only even indices. An exttrenum of 

F(I^,I^,*..) may be eimressed as 9F/9I « 0 for all n. 

o £. n 

Substituting (109) into (106) and using the ortbogonallity con- 
dition for the f^ (47), we obtain 


where 


Fd^,!^,...) - 

z I- I 

2 n n 


n 


fi 

N(I^, 120.0 ) •» Q 

S(x)(I 


** • ( 110 ) 
I^)[a^0(I-I^) -h ,(111) 



)i 







u 


It 1 » tmderittood that to b« eobstitutedi for I In ( 111 ). 

The condition that the 3 F/ 3 I^ vanish sieulteneoualy Isada to 


Vn “ Q / S<*) «' X)1<S* 


• Qh (k ) , 

n e' • 


( 112 ) 

(113) 


tfhere we have ua«d ( 29 ). Kote that ( 113 ) represent# the decos^ositicit 

of the energy balance equation ( 27 ), 

In the neighborhood of an extremna let the radiation field 

corresponding to a particular extretaum (steady state) be given by 

l^°^(x), or equivalently by I Nearby (in function 

space) we my write I(x) •» I^®^(x) +^(x), or the deviation ssay be 

written in terms oT its spectral components (> Expanding 

O 2 

F[X] about the local extremum, we obtain 


Fd^.I^....) 




Vm — 

(U4) 


the subscript o denotes evaluation at the extresam. The teras 
linear in 4 ^^ vanish bv.cauae 3 F/ 3 I^ ’^aniahes at the extremusa. Up to 
the tensis consiidered,^ F is locally a c|uadratic in ^ • The matrix 


tm 


(ai 91 I 

\ n m / « 


( 115 ) 


are the structute constants for the geometrical aurtace, .♦). 

If all eigenvalues of are positive, the surface la concave upwards. 
If one or ta^'re of the eigenvalues are negative, the surface la locally 
» atddle point, W« proceed to show that these eigenvalues ere the 
etabillty eigenvalues studied earlier. 





V* > - o ' s, . ’ 

^ ... 
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First note that if the tesperature field is alloited to b® « function 
of tits®, then by foil wing the approach for the 0-dln«Bion®l laodele we have 


X. 


fr 7 \ 

“’§r* (lo»^2****^ • 

n 


(116) 


i*e.* the tiM derivative is given by the gradient in this Bultidissen^ional 
sp«v«. For InfiniteBiMl departures item steady state m set lX„(t) " 1^** 
+ in (116) ( expand about ^ " 0» and obtain 


» \ n » /o 


tm 


(U7) 

(1X8) 


The latter equation concludes the proof that the local uoometrical atuctuta 
constants of at a steady state solution y>^ld the stability 

eigenvalues for that particular t^temdy state* Equation CdliS) Is the analog 
of the sliaple equation (80) » or in the one^dixsensionaJ equation (9T), 

Finally* as a conclusion to this section* consider the time behavior 
of the value of F(I^*l2t..*> vifteu the point (l^^^****^ ^ governed by 
the tltae-^depend^nt equation (89): 


dt 


■V ’>« ” ■ 


(119) 


( 120 ) 


wt»erc we have used <116). This latter 1« the multiditwnoional analog of 
(85). It has a cof<“«pondlng Interpretation: initial dbpartures of 

the state (1^,1^,.,.) frots a local extrensute of F lead » a trajectory 
of the point such that V decreases. The point will conrlnuo down the 
gradient of F until a local estrem?8 Is fotmd. Clearly local maxlssa 
®td saddle points are lassteble.: ; ; . : 

It was shown in Section 2 of this paper that a twe-eiode truncation 
of the apeetral representation gives a good approxlnaiEim in tieny 


ceeee. In this case vs cm actually plot a contour of F in the 
X ,1- pl(ssi«. This shwffl in Figure 1$; the exa^le is for constant 

O 2 

diffusion 80 that f <«) « P (a), Kote that the scale of and I» 
n n o i 

ere different In the test eioca the P (x) are not crthonoreal. This 

n 

two^di^asional imp :1s to b« with Its one»disa©ns tonal analog 

(gero-dissicivs tonal mod^l) in Figure 17* The labels I* II and III 
correspond to the analogous labeling of earlier figures* Physical 
interpretation of the figure follow® directly fro® the previous 
discussion* 

Before leaving this subject » we note that the Lyapunov functional 
Fix] gives rise to a variational principle for the Budyko-Sellers 
climate eodels* Although the functional developed here was strictly 
a laathematical construction* it Is tempting to speculate that there 
is an underlying physical principle analogous to extremum conditions 
in thenaDdynamics* If FlTj can be related to the rate at which 
entropy Is dissipated in the system* then the extremum principle 
would be in line with Frigogirte*® theory of nonequilibrium thermo- 
dynamic states (Frigogine* 1968)* Golitsyn and Hokhov (1978) have 
shown that the linear climte models (no ice-albedo feedback) can be 
formulated so that they satisfy Prigogine*8 condition* Unfortunately * 
such models have but one stable solution and* therefore, fail to test 
the general appUcabiUty of the extressum condition* Paltvidgc 
(1975, 1979) has suggested on the basis of energy budget observations 
that the climate governed by ® soiaevhat different* but 

equally Eimple, frinciple* h-hether the behavior of such a 

co?r 4 ?Iex can be bo ly characterized is a tantalizing but 

open quest ion* 
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7» Stationary perturbation theory 

The sensitivity of cllsats isodels to at&tionary perturbations has 
received consldeireble attention. Many recent papers have offered 
estisMtes of tha change In surface tet^erature caused by changes in 
the concentrations of varioua radlatively active trace gases in the 
atssosphere, by changes in th® surface albedo as a result of ttao’s 
activities, by changes in th® aerosol content of the atsiosphere. by 
additions of va»te heat from energy consusption, end taany others. 

In this section we shall derive a formula for the Infinitesimal 
change in climate as the result of small perturbations added to the 
energy balance equation. The formula has little practical utility; 
nevertheless, It offers insight into the qualitative response of a 
large class of models. For simplicity we shall study perturbations 
which are sufficiently staall that linear approxlsatlons can be 
applied. Analytic solutions to nonlinear systess with large per- 
turbations have been obtained by Salman {Saltaun, 1979; Salaua et al.. 
1980). 

Consider the class of energy-balance models discussed in Section 2 
There we derived a Green's function. given in (48), that describes 
the response to a ring of heat added at a given latitude, the ice 
line being held fixed. If G^ is known, the ice line latitude is 
determined for any Q through (50). In turn the temperature field 
given by (49) is cotspletely determined. Now we shall hold 0 • Q 

o 

fixed and add an ausount of heat g with some given latitudinal distri- 
bution u(x) (hetalsphevicslly ajiaBtetrlc, positive and normallted to 
unit integral). To add this heat (39) is replaced by ; ^ 
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L[T J + A(x) - QjjS(x) «<s,Sg) + g «(*) , 


(X21) 


end the integral equation analogous to (49) becos»s 


T (x) “ dy G^(x»y) S(y) e(y,x„) - A(y) + g u(y)l * (122) 


Setting T (x ) " T in (122) gives an expression which detenslnes the 

g ^ £1 

ice line latitude for a given g« By equating this expression to the 
previous T expresslont (49) evaluated at x • x » we find that the 

8 9 

terns involving A conceit giving 


(Sa)^ Q(*b> “ (Sa)g + (u)^ g(x^) . 


(123) 


where 


(Sa) - j dy G^(*g.y) S(y) a(y,x^) , 


(u)g • j dy u(y) • 


(125) 


Hiis result allows us to determine x (Q ,g) from the unperturbed result 

s o 

*s(Q.O). 

We now assume that g is snsall, so that the temperature field will 
be nearly equal to the g ■ 0 value. That is, 

T » T + 6T . (126) 

g o 

where is given by (49) with Q “ Qg oud d'f is siaall. Correspondingly, 
the Ice line latitude is given by 


X^ “ . 

o a 


( 127 ) 


vhe;« is th« wperturbed Ice line end Is estall* B>y expanding 

«(y»*g) in (122) and QCx^) in (123) to first order in 6a^«» we aay 

T ^rid Q to obtaiti 
o o 

6T(k) • I dy G^(a,y) [ S(y) 0(5x 

^ *o J * 


( 128 ) 


dx (u)^ « 

*-. ■ df ^ 

According to (129)* the shift In the ice line for a given s la directly 
proportional to the ice edge sensitivity dx^/dQ. By substJituting this 
shift into (128) we obtain the final result 


6T(x) « g dz C(x,z) u(z) 4 0(g"^) * 


%diere 


r,(x.d) = 0„C..r)+^ f dy C„(,.,) S(y) ^ (y..,) C„(. ,.) 

U o o 


Setting u(z) « 6(z - zji in (130) shows that G(x*r^) represents 
the temperature response to a ring of heat added at a glvem latitude. 
The response includes to first order the effect of the Ice line 
shift. Since the problem has been linearized* the responses to an 
arbitrary distribution is given by the appropriate superpos?.ltion of 
localized sources* as is indicated by (130). 

Ve Should recognize that the first term in (131) is thte response 
when there is no ice-albedo feedback. The second terra is t;he responoe 
due to the albedo feedback* Kota that the fi&edback term is3 multiplied 
by dK /dQ* the slope- of the ice line aolar constant curve. TIuia* for • 
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sisall perturbations, the slope of the ice line curve governs the theroal 
response of a cllmte nodel to any perturbation. As a result the 
fundamental sensitivity (J^ - Q/iOO dT^/dQ of a »odel la an Indlcatox 
of it® sensitivity to my perturbation* 

For the step function albedo 6^ and dx^/4Q Esay be related by 
differentiating (31): 


io5-» * * »To ♦ 0^ S(.„) ^ 


(132) 


where Aa is the change in albedo at x " x^» (This la a particular case 
of (105)). As a mjaarical exas 5 )le, consider the caae of a etep function 
albedo for which the expression for G simplifies to 

G(x,z) - G^(x,z) + - °o.42‘C ®o^**V • (^33) 


where (132) has been used. (Recall from (9) that 1.12*C is the sensi- 
tivity when the albedo remains constant.) Since G^(x,z) peaks at 
X • 2 (see Figure 7, for example), the feedback terta tends to increase 
the response near the ice line, x == x^, and this effect Is largest 
when the added heat is closest to the ice line, * - x^. Figure 20 
Illustrates these features for a tropical source (z - 0,4) and a 
midlatitude source (z " 0.7), where we have used the diffusive G 

o* 

and x^ ■ 0.9, For a source which is broadly distributed 
only the peak at the ice edge appears. This kind of effect has also 
been observed In detailed models having many feedbacks (Hanabe and 
Wetherald, 1975), 
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Floes of g(x ) show that the rraige of g for which tlhe linear 
s 

epproxlTOtion Is valid tends to be smaller when heat is added closer 
CO X • This range also depends upon the Q(s ) curve characteristics 

O 0 

(SaliBun, 1979; Saltdun et al« » 1980). For a range of large negative g 
the slope, dx^/dg. usually changes sign and it can be sb^m by arguments 
analogous to chose in Section 6 that a negative slope is necessary and 
sufficient for instability. 



09 Fluctuations 



In previous auctions w have focused <m quasi-pariodic changes 
associated ?d.th seasonal and orbital However* typical cliioatic 

t±vm series also esdiibit large amounts of TOcperiodic variebility* 

In this section we shall e:stend the models considered so far by 
incltsding this broad distribution of variance having periods that 
range fro® seasons to centuries « The low frequency pairt of this 
range, say periods from decades to centuries, is particularly 
difficult to observe because instrumental records tend to be too 
short, and proxy records are difficult to interpret* Nevertheless, 
such low frequency natural variability my contain information on the 
true sensitivity of the climate, Toat Is, the effects of possible 
future perturbations of the clitaate due to \^rlous extern?.al causes 
may be estimated from past responses to natural fluctuatf^ons associated 
internal degrees of freedom. In Sections 2 and 6 thm stability 
of a given steady state was related to the sensitivity. One of the 
goals in this section Is to show hcr>g the Ut;tural variabllf.ty may be 
related to the sensitivity. Such a relatiorishlp was sagg-ested by 
Leith (1975, 1978) and is kncr«m in statistical mechanics ms the 
’’fluctuation dissipation theorem.” 

In long-term climatic records slower fluctuations gernerally have 
larger amplitudes. This feature is also a characteristic vof Brownian 
taotion, the random moves^nt of small particles suspended im a liquid. 
Over short time periods only a few saolecules collide with a given 
particle and Its displacement frost its original location ils small. 



but over long periods there Is eoffi® chance that many collisions drive 
the particle far from its original position. The fluid has two 
coBjponentss the selecules, which aove and change direction rapidly; 
and the particles j, which laove slowly. A sialiar scale separation 
Is expected to exist between the relatively rapid evolution of 
synoptic weather systems and the more sluggish climatic cotaponents 
such as the global average temperature and the pole-to-equator 
temperature gradient. 

The suggestion that broad~band climatic fluctuations may be a 
cumulative effect due to variations on much shorter time scales was 
made by Mitchell (1966) in connection with sea-surface temperature 
anomalies. The two-tlme-scale approximation was given a general 
formulation by Hasselmann (1976) > who eim>haslzed the role of negative 
feedback processes In limiting climatic variability. This approach 
was applied to simple ocean modsls by Franklgnoul and Hasselmann 
(1977) and Franklgnoul (1979), to Budyko’s energy-balance model by 
Lemke (1977), and to a global energy-balance model by Fraedrich 
(1978). Sobock (1978) has performed numerical computations of 
effects of fluctuations in energy-balance asodels. 

a. Preliminaries 

Me shall think of the coc 5 )lete time history of the global 
average temperature, for e>;ample, as a single realization of a random 
phenomenon. If we imagine an infinite population of essentially 
identical planets in the saiswi orbit around the sun, the collection or 
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enseTTible of all the time series of global teasperatures constitutes a 
Stochastic process, generally defined as any ordered (by ti®e in this 
case) set of random variables (Jenkins and Watts, 1968). The proba- 
bility that the temperature is in a given range may be estimated from 
a saisgjle of the population by determining the fraction of planets within 
the sample having toaparatures in that range. Such a result generally 
depends on when the count is taken, particularly if there are changes 
in external forcing. We shall assume that such dependence aay be 
removed by subtracting the effects of external forcing. In our 
example of an ensemble of planets, this might be done by considering 
only the deviations from a time-dependent ensemble average temperature. 
The associated stochastic process is said to be ,stattpna_^. 

Since no such ensemble of planets is available, we are forced to 
try to determine the statistics of any stochastic compcnents of the 
climate from a single realization. This might be possible tf the 
values of climatic quantities at a given time are uncorrelatad 
w-:ch the values at a much later time, so that the complete record may 
be treated as an ensemble of independent records. We shall assume 
that this is the case, and that any quantity determined by time 
averaging over a sufficient length of a single record will equal its 
corresponding ensemble average value. Ihe stationary stochastic 

process is then saici to be ergodic a 

A general stochastic process, then, is an ensemble of ifcictions 
that depend on time. Each function represents a point in a "sample 
space" of possible experimental outcomes (different planets in the 



above esasEple) which we ®ay label by the paraetster e« We use the 

notation T^(t,e) to represent a different function of tlsws for each 

value of e. Note that T (t,e) is also a randoa variable in e at any 

0 

given tisae. We will use angular brackets to denote an average over 
the enseable. The tasaa, for exaasple, is given by 


<T <t)> H lira ~ 


N 


i T„(t,e) , 


j5^ ” e«»l 
and the eutocovariance function is given by 


( 134 ) 


<6T^(t) 6T (t + d)> 
o o 


“ lim 

JJ-WO 


1 

K 



6T^(t,e)6T^(t + T,c) , 

( 135 ) 


wh'iJ'e 6T Is the deviation from the mean. Presuming the effects of 
o 

external forcing have been subtracted^ these quantities will be 
8tationar>s i#e*. Independent of t, 

Tliere are two related techniques for detenaining the statistical 
properties of a stochastic process, namely the Fokker-Planck and the 
Langcvin methods * The Fokker-Planck method deals directly with the 
probabilities of various events, and characterizes the type of process 
through relations between these r tobabilitles. The Langevin method 
begins with a deterministic equation for the mean motion, adds a 
atochastic forcinR with assumed statistical properties, and uses the 
resulting stochastic equation to derive the statistical properties of 
the motiotu We shall employ the Langevin method sxnce it requires 
only a simple extension of our models for the mean clltsate*. 



5.^i5g lg_modal with stochastic forcing 
To the extent that the various nodal parameters fluctuate on 
time scales typical of ••weather,” l.e., well-separated from climatic 
response times, we will assume that their collective effect is that 
of a random forcing. Of course, the statistics of such forcing 
should be determined from detailed models such as general circulation 
models (GCMs) or from meteorological data. As a first guess, however, 
we will take the amplitudes of the forcing to be distributed according 
to a Gaussian distribution* 

The stochastic version of the simple global model discussed in 
Sections 2 and 6 has the form 
dT 

C ^ + A + BT^ . Q H^(T^) ^ 

where 6n represents the random forcing, measured relative to the mean 
forcing, so that <6H> « 0. The expression (13 d) represents a different 
equation for each member of the ensemble, l.e., for each e. The 
ensemble average of all these equations gives an expression for 
<T^>, and we shall assume that <T^> is one of the stable steady 
states considered in preceding sections. 

If the fluctuations in T^ are sufficiently small, we may linearize 
(136) about a steady state. Upon linearizing (136) we obtain 

^ A . ^T4*.x<st 

dt o 00 oovt.t; , (137) 

where .= - <T^>, and is given by 
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with fixed albedo 
with albedo feedback 


( 138 ) 


where S Q dT^/dQ Is the sensitivity. According to Section 6, \ 
taay be Interpreted as the stability eigenvalue of the mean climate 
<Vo>. It is given by the curvature at the associated minimum of the 
global model potential shown in Figure 17. The formal solution to 
(137) is given by 


<STo(t,e) - f^dt’ g^(t - t') 6H(t’,c) , 

where 

8j,(T) « i 0(T) e ° 


(139) 


(140) 


is the retarded Green's function for the linearized equation. The 
unit step f, action 0 (t) in (140) vanishes when T < Q. As a result 
the response at time t depends only on the forcing at earlier times 
t’ < t. 


Note that the linearization requires only that remain small. 
The forcing can be quite general. Thus, we may apply (139J to 
determine the response to deterministic changes in external heat 
sources, as well as the response to fluctuating intenial beat sources 
For example, substituting a Dirac delta function for 5H shows that 
g^(T) may be interpreted as the effect of a heat impulse at time t' 
on the temperature at time t' + T. For a constant change in heating 
beginning at t’ >» O, the temperature response at time t is the 
integral over g^(t) for t ^ t. This integral grows to a final value 
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proporf tonal to 8^ after a chf .’acterietlc response tine given 

by 


C « ^ ft 
^ » ■ 


(Ul> 


So the response ti«» m «ell as the shift Jn the taean temperature are 

both detersilned by the sensitivity. 

The KtatlHticaX properties of randoo variables Bay be determined 

either by specifying their £roba]^t^ specifying 

their moments (see. for example, PapouUfi. 1965. Section 5-4). If fiH 
mmojsa a Cau-sslan distribution, its odd moffiopts will be zero. 

and Its even ca^taents will be given by a sua of second moment products. 
<6H(t.) 6ti(t-,)...fiH(t, )> " ^ <6H(tj) 6ll(t.)> <6H(tj^) 

I d. *•» f.\% 


ell 

pairs 


( 143 ) 


As a 


result, the atatietlcs of «H are completely described by the second 


moment, cr autocovarianct. which we denote by We shall assume Chat 

Che autocorrelation tlsas of the forcing «ay v.e neglected relative to 
the response time of the temperature so that we may write 


r*’<T) 5 <5H(t) fiH(t * T)> » y 6(f) , 


( 144 ) 


where T is a constant. Tl,i. Is called whUe nolBo in analogy to light 

i,\ € coTtx*Bi*onCAnz ?jpeccrvm* \dhleh is given by 


Ji 




dt T^*(t) * Y , 


( 145 ) 


is i«dependi?nt of «, the angul-;r frequency. 
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Using the linear escpresslon (139) and re®s 9 ©barlng thast can be 
factored from any enaec^le average, one can ehov that the' 
of 6T also satisfy (142) and (143). Ihus, 6X is also Gmselan, 

O 0 

and its acaeistica are described by its autocovariance, uhtlch we 

m m 

denote by , Applying (139) gives as a double Inregral 
over . Using (144) we obtain 

{t(/C 

r*„(T) H <6T (t) 6T (t + T)> - (Y/2X C) e ® , (146) 

w V iftf 9 w 

Kote that the characteristic decay tl<s« of the autocovarlarace, which 
is the autocorrelation time Is identical to the reoponste titw for 
a step chaTiga in forcing , 


C/K 


(147) 


The autocorrelation time is thut also d<itemlned by the oeiasltlvlty. 
The relation (147) ia a special ease of the tsore general fluctuation- 
dissipation relation given by 


(T)/r'^ (o) • C g (t) , T > 0 . 

oo oo o 


(148) 


tlie left sid§ of (148) is potentially measurable from cllBS&tic records* 
Tnrough (139) tlie resulting knowledge of vouXd allow the estimation 
of global, tesaperature v‘hangea tha; might result from any ccff/jtetr?»lated 
change i... global heating provided we could estimate the thex‘faal Inertia 
C* l>r a step change ±ix heating we have seen that the res$H>nse 

approacheis the integral of g over sll T and^ Binc^ both of 

. ■ ^ 

(HB) ijire cxpo^tentl^5la in T*, the lnt<^grmX. of, (148) reduce© to (14?)* 
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We refer to (148) «s the ''Rtrong*' fora of the nwctuaeion-^laslpetlon 
relation. It detenainee the response at any tiiae T toj any given 
forcing. We refer to the integrated version [equality? of the integral 
of left and right-hand sides of (148)1 as the ’Wak" fcor®. It deter- 
mines the response at large time® to a step function itn forcing. The 
weak fora laay retsvaln a good approxissitlon In forced diisslpatlve systems 
even If tVie strong fora does not hold at all T (Bell. j1979). 

Figure 21 providef. an Illustration of the flue tuattlon-dlsslpat ion 
relation. Shown in Figure 21a is a typical time setiem of T^. The 
time series could be either measured or modcl-generatedS. It contains 
both spontaneous fluctuations and a shift in the mean ast t » 0. 

Figure 21b shows the taean tenq^ernture. The taean may bt© extracted 
by time avei aging segments of the record shown in Figurre 21a. If 
the record were model-generated, the mean would be obtaalned by 
averaging an ensemble of model runs. In time the mcean temperature 
responds to a change In mean heating. Subtracting Figt^sre 21b from 
21a prodisces a stationary time series for 6T^ “ ~ *^^To^ 

zero mean and an autocovariance with exponential form. ms is shown in 


Figure 21c. As suggested by the f luctuatlon-disslpatlo?n relation, 
the autocorrelation time (indicated In the figure by Tf^) is identical 

to response 

The above results imiy be re-expressed in the freqiaency vJnraain by 
Fourier transforming. Transforming (140) gives 


o 


1 


iiX X. 



# ( 149 ) 
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and tranafoming (146) givea th© variance epectrtiss for the global 
tei^erature 

sjp(w) “ |g^,(«))|^Y - ■ TY ' ^ "" r2 ' • (150) 

” W^C + ^ 

O , 

Again 0 ) Is the angular frequency* This epectrua Is ache&atically 
illustrated in Figure 22* *^:;he figure exhibits a feature cecsBon to 

aany clitaatlc tliae seri'^a* nesely the variance increases as the 
frequency docreases. This concentration of variance at low fre- 
quencies is tenaed "red noise." According to the linearized feed- 
back ffiodela, the spectrim grows as uT^ and flattens below a frequency 
proportional ^ that this frequency is inversely 

proportional to the sensitivity* 

The tmcertainties in climate sensitivity discussed in Section 4 
imply uncertainties in the autocorrelation time of the global tem- 
perature and in the corresponding frequency dependence of the tem- 
perature spectrum« According to (138) and (147) > T would be equal 

c 

to the radiative relaxation time C/B If albedo feedback were negligible 
If this were the case, the spectrum of the global temperature would 
Increase as <4T*“ <jcywn to frequencies near B/C~s few tenths of a cycle 
per year. On the other hand, if albedo feedback is large, the 
stability parameter becorses small. If the present climate is only 
marginally stabw^e (such as the state labeled I in Figure I, corres- 
ponding to the shallow minimum in Figure 17), then the outocorrelatlon 
tlm Bay be csich longer then C/B and the w *” behavior may persist to ■ 
smeh lover frequencies. According to (150) our ability to reduce 
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uncertainties in clltaate sensitivity vflll depend on our ability to 
estimate cliaace spectra at low frequencies. 


c. Zonal cllmato tRodels with, stochastic fgr£i jl& 

As for the sisple global mean saodel, fluctuation di-sslpatlon 
relations can also be derived for tonal laodela. In tonal nodela» we 
expect the response to fall off over aone characteristic distance, as 
illustrated for the diffusive oodci in Figure 20. According to the 
fluctuation-dissipation relation, this distance npproxloately equals the 
correlation length derived from the spatial cross covariance. Regions 
aeparated by ®ore than this characteristic distance represent lndepende.:t 
climates. As a result we could, for example, apply the same global-type 

seodel to two such independent regions* 

We shall restrict outselv* > here to tc>nal average Ts^ dels with 
constant C* To study the zonal average models we will follow the 
approach of Section 6* For a basis set we will use the eigenfunctions 
of the linearized energy balance equation* With this ck-pansion, the 
potential cur\»ature of the global model, will be replaced by the 
eigenvalues of the linearized equation, As in Section 6, the 

are the curvature parat^ters of a generalized potential* Tlie autCM- 
correlation times for each ® decrease as n iticreasea* 

Consider the class of onc^dimenslonal stochastic models obtained 
t ' f-iding random forcing to. the tit^^dependent , version of (39) r . 
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C + LtTl » H + «H . (151) 

As before, w« asstasa 6H to have zero asan end“ to be uotiaally distributed. 
We wlU also take the spatial cross-covarlanco to be given by 

<(SH(x,t) 6H(y,c + T)> “ Y(x,y) 6(t) , (152) 


where 69 before the brackets indicate an ensenble average. 

Linearizing (151) about a atable steady state gives 

C |- 4T + Lt6T] - H’16T) + 6H . (153) 

where the term H* Is evaluated for the state In question. 

In order to relate the fluctuations to the sensitivity as we 
did for the global model, we assume that L - H’ has a discrete spectrum 
as pictured in Figure 18, and expand 6T In terms of the corres- 
ponding eigenfunctions ss follows: 


fiXCx.t) • 

u 

With a similar expansicn for 4 h, the Individual modes will fluctuate 
according to 


C 


dt 


JT + X 6T 
n n n 


6H 


n 


(155) 


and these ere coupled through 


<5Hp(t) (t + l)> “ 5 (t) 


(156) 


where 
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fl rl 


0 'o 


dy ^^(y) y(n,y) . 


Since (155) is analogous to the global form (137), the autocovariance 
and spectrum of the are analogous to (146) and (150): 


r^(T) 5 <6T^(t) 6T^(t 4- T): 


•A_t/C 


X T/C 


6(T)+e “ e(-T)l .(158) 


S^(u)) 


(IwC + X ) (-iu)C + X ) * 

ui n 


Modes with successively higher values of X^ (smaller spatial scales) 

hove less variance, shorter autocorrelation times and spectra which 

tend to be flat at frequencies such that uT^ 

oo 

As with Che simple global model, (153) may be solved through the 
use of a retarded Green's function. The Green's function satisfies 

(c|^+L-H*) G(x,x'; t-c') - 6(x-x’) 6(t~t') , (160) 


In terras of the Green's function the temperature fluctuatl 


ons are 


given by 


6T(x,t,e) • / dx' /dt' G(x,x'; t-t') 6H(x',t*,e) , (161) 

By using (161) and (152) to evaluate the spatial cross-covariance, and 
by performing one of the tltae Integrals, we obtain 


r <x,y;T> *. /dx’ /dy' /dt’ G(x,x';t') G(y,y';t' < T) Y(s’,y') . 


\ . 
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This fomal result is independent of the basis functions. Note that 
th® Biodal expression given by (158) nay be recovered froa (162) by 


setting 




where 


-X T/C 

g^(T) “ ^ e 


From the orthonormality of the it can be shown that 

G(y,y'{ t’ + t) ■ C / dz G(z,y*jt’) G(z,y; t) • 


(163) 


(164) 


As a result, (162) becomes 


r'^(x,y;T) “ C / dz F^(x,z;0) G(z,y;T) • 


(165) 


T “1 

If we define the Inverse of the zero-lag cross-covariance, [T ) , 


such that /dzir^r^(x.z;0)r‘'^(t,y;0) - «(z-y). then (165) gives 


/dz[r'^r^(x,z;0)r‘^(z,y;T) - CG(z,y;t). 


(166) 


The expression (166) is a generalization of the global fluctuation- 

dissipation relation (148) and it is the ’’strong” fora since it applies 

for each T. It gives a direct connection between the unpredictable 

T 

natural variability of the steady-state climate, as measured by T , and 
the potentially predictable response of the mean climate to an 
impulsive external perturbation. According to this ’'fluctuation- 
dissipation relation," if the thermal inertia of the system C Is known, 
an estlisace of from the climatic record determines G. Using G, 
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the 8»an response to any perturbation, the addition of vaste heat, 
a change In albedo, etc*, is obtained at any given latitude and tisKs 
by convolutlng the proposed change with G, as in (161)* Consider, 
for example, a constant change In heating added at x* y and beginning 
at t' •• 0. From (161) this change produces a response at x and at 
time t given by the Integral over G(x,y;x) for 0 ^ T ^ t. For large t 
this integral approaches the stationary response function, G(x,y) 
discussed in Section 7, and Illustrated for a particular model in 
Figure 20. According to (166), the stationary response of the true 
climate to a perturbation could be estimated from the Integral of 
r (Xty;T) over all positive T. As for the global laodelg low frequencies 
dominate, and in terms of the spatial cross spectnmi we obtain 


/ dz[r’‘]"^ (x,z;0) ~ S^(z,y;0) + i f*~ 

, c 


C G(x,y) 


(167) 


In our terminology (167) is referred to as the ”wcak” form of the zonal 
fluctuation-dissipation relation. It generalizes the global version 
(147). The zonal version not only relates the time scales of natural 
fluctuations and responses to perturbations, but their spatial dlstrl- 
butiOins as well. According to (167), the s\*>atial response may be 
estimated from the low frequency cross-covarlanr e divided by the total 
cross-covariancc. 

In closing this section we should recall our original assiaaptlon 
that the system can be divided Into slcjw and fast components, i.e., 
that the autocorrelation time of ‘Wather'V is mich shorter than climatic 
.s. Justification for Mitchell’s original stochastic 
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treataant of n&a eurface temperature anomalies Is provided by the 
relatively long radiative relanation time of an ocean coltmsn compared 
to that of the atG»:spheric colum that provides the forcing. In 
8udyko~Seller8 models the ice~albedo feedback results in long response 
times at least for the largest spatial-scales, but one must also 
isolate the effects of the mean forcing in each mode. The stochastic 
treatment is less tenable for the higher modes since, as we have seen, 
the response times tend to be shorter and the mean forcing less 
well-understood. 
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9» Disct;tai3ion 

It is remarkable that the tonally averaged eea level temperature 
for the northern hemisphere can be fitted to the simple parabolic 
form 

T(x,t) ® cos(2Trt - 4».j.)Pj(x) + T 2 ? 2 (*) (168) 

with an rms error of only 2*C« The parameters T^» Tj^, and T 2 thus 
give a reasonable representation of this hemispheric climate variable. 

An obvious goal In global climate theory Is to construct models that 
generate these four parameters In agreement with observations and pre“ 
diet how the parameters change when external conditions change. 

Since (168) involves only the largest space scales and the annual 
time scale and since the solar heating Is also describable by these 
same scales, we are encouraged to try simple heat balance models that 
connect the two. This review has surveyed and appraised recent progress 
toward achieving the connection. 

In Section 2 a sequence of mean annual models was introduced. 

It was shown that the mean annual global temperature T^ can be esti- 
mated using the most elementary radiation balance considerations 
provided an empirical formula for the terrestrial outgoing radiation 
is used. This latter includes corrections for the greenhouse effect 
due to the presence of infrared absorbing gases, clouds and even the 
change of absolute humidity with t< erature. Provided the observed 
present albedo is used, this e<!tlmate of Is independent of trans- 
port mechanisms. The sensitivity of the climate 6^ was shown to 
be about twice large for the empirical earth model as for an 



ideal infrared ealtting planet. The enhancement can be accounted for 
by the change in absolute humidity with temperature. 

In order to estimate T^, which is a measure of the mean annual 
pole-to-equator temperature gradient, one wist introduce assumptions 
about how heat is redistributed on the earth’s surface by the geo- 
physical fluid system. Because diffusive terms lead naturally to the 
parabolic component P 2 <x) in the energy balance equation. It is natural 
to take the transport to be proportional to the gradient of temper- 
ature. For simple diffusive models, the diffusion coefficient may 
be adjusted so that the correct amplitude is obtained. In order 
to estimate the sensitivity of the resulting two-mode model, more 
assumptions must be Introduced. For example, if the earth is cooled, 
the ice caps expand, the planet becomes more reflective to solar 
radiation. To allow for such changes, feedbacks such as cloudiness 
changes, changes In the transport model and changes in the Infrared 
formula due, for example, to lapse rate changes ought to be Included. 

As such feedbacks remain largely unknown, they were not created 
extensively in this paper. They, of course, merit further attention 
as more empirical and theoretical Infonaatlon accrue. 

In Section 2 it was shown vhat the simple ice-cap models are non- 
linear and have solutions that exhibit a rich structure. In particular 
there are at least three solutions for the present value of the solar 
constant. If the sun’s luminosity is lowered by about 10%, the model 
climates experience a catastrophic transition to an ice-cov«red earth. 
This feature appears to exist even in the results produced with 
nufljerical general circulation cllia-ate models. ' ■ . ' r; ■v;? 
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By introducing a heat storage term C3T/3t to the energy balance 
eqimtion and alloying the solar absorption to vary with the seasonal 
cycle, we showed in Section 3 that and could be estimated. It 
appears that reasonable values for these parameters can only be obtained 
if the value of the thermal inertia C is taken to be much staaller over 
land than over ocean# As a test of the single seasonal model con** 
struction, the four parameters T^, and were computed for the 
southern hemisphere by changing only the albedo and the land fraction. 
The test proved satisfactory enough to suggest that the main feed- 
backs which operate on a seasonal time scale had been included. 

Section 4 concentrated upon the relationship between tbe quality 
of the parameterization f rmulae and their effect upon estlmtes of 
the sensitivity to solar constant changes. Each empirical formula 
was criticized and the need for further work in parameterlzations was 
emphasized. 

Even with the reservations of Section 4, we think that some 
problems of practical interest can be examined with the sixaple models. 

In Section 5 we used the seasonal model with ice-cap feedback to 
estimate the effects of changes in the earth’s orbital elements. We 
found that the climate response to orbital changes was about an order 
of magnitude less than the paleocllmatlc data seem to indicaite. We 
suspect that this is a large enough discrepancy to rule out the sitnple 
instantaneously responding ice-cap feedback as a prime cause of the 
ice ages. We emphasize that many lower frequency feedbacks cannot be 
discover<^d by. , examining only the seasonal cycle. Clearly, , trore work , 
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can be done on this problem and the simple models will provit’e a yssiul 
framework for future discussions. 

Section covers the stability of the simple models. In parti- 
cular * it was shovm that the linear stability of a solution can be 
deduced from the :lgn of the local slope of the ice line versus solar 
constant curve. This slope-stability theorem appears to hold for a 
broad class of models Independent of numerical inputs. An analysis 
is also given for the finite amplitude time behavior of a class of 
models, A potential function can be constructed that yields quali- 
tative inforniation about the behavior of the system far from steady 
states. Tlie potential function can also be thought of as a vari- 
ational principle fcr the climate. 

Section 7 presents a simple analysis of sca'^lonary perturbations 
to the heat balanc^-. The analysis shows that the sensitivity for any 
perturbation is related to the sensitivity for solar constant c''"angcs. 
The relationship is shown to hold for a broad class of models. 

Section 8 concludes the review with an Introduction to stochastic 
clima'ui models. The sitr.>le models of previous sections were linearized 
and allm<fed to have a stochastic white noise forcing. Analytical 
solutions are easily found for the various climate statistics. Example 
proofs of the fluctuation-dissipation theorem were given. Through 
this theorem, it was shown that Inforc-ation about climate sensitivity 
can be derived from data on natural fluctuations. 

Hiiny of the results of this review are expected to hold for more 
cotiiprehfinsive climate models. Hie simple models foim an intuitive 




94 






froa vhich ws ean study the larger but less scrutable eiodela* 

For eKsapIe, vhen a result fails to carry over to a larger model* 

%m Isaadlately face the problem of discovering «hy* Future experi- 
sauits with the large and small models my lead to valid perameterl- 
sations that can be used in the simple models* If this turns out to 
be true* the simple isodels my be able to cosdbine the psre^teriaation 
formulae found vith empirical date in such a vay as to pley a sig> 
tuficant role in the aasessiasnt of climte change- 
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IkPFSKDIX X: Hotfttlon and Hmserical Values of PeratMters Used 
Equation nus&ers given refer to the first use of the syisbol 
in an equation or to e^iuations nearest the first occurrence of the 
eytdtol in the tent* In some eases a nuaerleal valw* is listed which 
coincides with that used in tha text* These nuahers are likely to 
change with it^roved taeasure&eats and laodele* 
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V/ 


••2 

A intarcept in the Budyko rodlation fozmla (7)* 203*3 a 

ioT nottharn heaiRphere. 

aCa) albedo at latitude designated by x (2)* Cf*» Table 1 for 
northern hemisphere* 

C(^ planetary albedo (2)* 0*30* 

^l*^f absorption fraction (co>albedo) of earchoattaosphere system 
over ice and ice*free surfaces (13). 0.38, 0.70, 
average co-albedo at ice-cap edge (18). 0.54. 

%**2 l<agendre coefficients of co-albedo (19). 0,681, -0,202. 

a(x,x^) co-albedo as function of latitude and ice-cap edge, 
a^^Ct) time-dependent coefficient of Pj(x) (61). Cf. Table 1. 

B linear coefficient in Budyko radlftlon formula (7). 

2.09 W n^^rc)"^. 

8^ sensitivity (5), Value tsodel dependent, as discussed 

in Section 4* 

C effective heat capacity per unit area of earth-atmosphere 

system (78), 

^ ideal land and oceanic nixed layer. 0.16 B yrs., 

4.7 B yrs, 

D thermal diffusion coefficient (22), (33). 0.649 W n“^(*C)'’^ 

or D/B - 0.310. 

affective D for latitude dependent diffusion In two-mode 
approximation (36), 

&a change in coalbedo at the ice cap edge (99b), 

e eccentricity of earth’s ovblt (73), 0.017 (present value). 
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f (k) orthonomsl elgenfuncClona of the dlffuelcn operator (47)* 

f, ,f fraction of land and water In a latitude belt fox the 

la W 

eloplified aeaaonal ^del (54) « 

F(T^)i Potential function for the sero and one**dimen8ional taodela 
(86), (1)5). 

FtT) Lyapunov functional for the tcnperature field (107). 

infinitesimal departure from a local extremun in F(Iq,I 2 «**) 
(115). 

S global average value of the stationary heat perturbation 

( 121 ). 

g (T), tetaperature impulse'-response functions for t;he global 
o 

g (T) average or higher modes (139), (163). 

D 

Gr«en*s function for the linear transport operator In the 
energy balance equation (42) ♦ 

G(k>ji:) Green’s function for the transport cosibined with linearized 
ice-albedo fee;Iback* 

G^^^ ©ode coupling coefficients (64)* 

G(xpX ;t-t’)C#reen’a function for the linearized tls^-dependent energy- 
o 

balance equation. 

Budyko transport coefficient (34); also the strength of 
global white noise forcing in (144). 

Y(x»y) spatial cross-covariance of the white noise forcing in 


( 152 ) 
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Y espgneion coefficients of Y<*»y) t*'® 

r*^<T) autocovasrlance of globsl white noise forcing at lag T» 

(■?) autocovarinnce of global tciss>erature fluctuations at leg 

oo 

tiee T« 

r"^(s 5 ,y»T) spatial croes-covariance of sonel tes;>erature fluctuations 
at lag tis3 T. 

(t) expansion coefficients of F^(x,y;T) in the basis. 

wx 

H (* ) Legendre polynomial amplitude of solar heat ebsorbed (29), 
Q 9 

h (^ ) Q^lityds of solar heat absorbed with respect to the 
n s 

eigenfunctions f (jt) (91b), 
n 

X,X(x)» infrared radiation to space (W «”^), (7), also Table 1, 
X(x,t) 


X 

8 

X 

n 


L ,t 
n* n 


X 

V 




P_(x) 

n 


1 evaluated at x ■ (90). 

Legendre or f^(«) amplitude for I (90). 
a linear operator (40). 

eigenvalues of the linear operator in an energy balance 
equation (29), (49), 
stability eigenvalue (98). 

an empirical transfer coefficient coupling land and water 
areas in the seasonal model (56), 

orthcnorcml eigenfunctions of the linearized steady^state 

energy-balance equation, with eigenvalues given by the 

stability parameter X. In the absence of ice-albedo 

feedback, <j< ■> f and X « i . 

n n n n 

Legendre polynomial. 
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Q»Q^ solar constant divided by four (8). Subscript denotes 
present value (335 W aT^ in Section 2, 340 W sT^ In 
Section 3). 

R radius of earth (1) ♦ 

O Stefan-Boltzisann constant (1) 0*5669 x 10*“^ W 

o solar constant, see Q,Q (1), 

o o 

S(K)f Theoretical distribution of solar host energy reaching 
S(x,t) top of the attsjosphere (2), (3), Table 1, (74), 

S^(t) Legendre aode aoplitude for S(x,t) (3), (61J, Table 1, 


(74). 


sLm 

iSn 


T, T(x), 


T 


R 


T 


n 


T 

s 




u(x) 


spectrua of random forcing. 

spectrum of global temperature fluctuations., 

cross spectrum of temperature fluctuations lisn modes m 

and n. 

Zonally averaged 1000 mb level temperature fitLeld. Wlien t 

does not appear, annual averaging Is Implied^ Table 1. 

planetary radiative temperature (1). 

Legendre i^de ao^lltude for T(x) or T(x,t). T is the 

o 

planetary average tetoperature. 

temperature at ice-cap edge (16). Budyko’s a-iile; 

- -10*C. 

lag time, climatic autocorrelation time, clteatic response 

time (values depend on 8 ). 

o 

nonsalired distribution of added heat flux (E21). 
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aim of Istituda (2). 

sine of latitude «t iceo^cap edge (13)* For the present 

cliitate» n • 0*95* 
s 
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APPENDIX II: Closed Fora Green's F-^ction 

In this appendiK we present a brief derivation or the Green's 
function for a diffusive cllisate eiodel* Instead of the fora (2«45)« 
which converges slowly* we seek a closed form espresslon. We follow 
a notation used by North (1975a) in a slightly different context. 
Consider the differential equation 

"I ^ ^ G^(x.x’) + G^(x.x') - «(x - x') (A.l) 

and subject to the usual vanishing gradient boundary condition at the 
pole and equator. For x > x* the equation Is hoaogeneous and is known 
to have the solution (Kamket 1959) 

Gp(x,x’) - A^P^(x) + A^Q^Cx) . X > x' (4.2) 

where P^(x) and Q^(x) are the Legendre functions and 

^ . 1 + i (1 - ^) 1/2 (A.3) 

V 2^2^'^ D' * 

which in general »ay be complex. The coefficient A^ oust vanish 
since Q^(x) diverges at the pole. Similarly, below the heat source, 

X < x' we find (Kamkc, 1959) 


Cg(x,x') 


A, fj,/*) + A^f2y(x') 


(Aa A) 


The functions P^(x), fj^(x) and f 2 ^(x) may all be related to hyper- 
geometric functions (Kanke, 1959; Erdelyl, 1953): 
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P^(x) 

. F(i + 1 V, ~ I V, 1, 1 - 



n y(«» i Vb 4* ~ V i 

(A.5) 

f2v<«> 

- xP(f - f V, 1 + i V. xh . 



Power series representations asay be used to evaluate any of «-hese 
functions in the doasain of interest. If ve make use of the property 

F( 2 ,b,c,z) ■ ~ P (a •{ 1, b + 1, c + 1, a) , (A.6) 

and 

^(®*h.c»o) “ 1 , (A. 7) 

it can be shown that imist vanish to satisfy the equatorial boundary 
conditione Now continuity requires 

Cp(x’.,-) - c^w.x') (*.g, 

which leads to 


A3 - P^(x»)/f^^(x') 


(A.9) 


Now by integrating (A.l) over an inflniteeimal interval about x - x* 
we obtain another condition: 


B df G(x.x') 


3c4- 


(A. 10) 


I 3C- 

which leads to the unique solution for Aj^, A^. The Green’s function 


may now b® written 
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(A. 11) 


(A.12) 


where the argument x' and the Indices v,l have been svE|»pressed and 
the prime on P and £ indicates derivative* 
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PIGOBB CAPTIONS 

Figure 1 Solid line depicts the fraction of solar radiation absorbed, 

Kq. (13), for the rero-diiBensional global clijaate 
taodel with variable ice cap. The dashed line shows the 
outgoing infrared energy per unit area per unit tlae 
divided by the present solar constant (H). Intersections 
of these curves arc roots corresponding to steady state 
cllmtes: Koot I is th*r. present cliiaate; Root II is a 
large lce~cap solution; Root III is a totally ice-covered 
planet. 

Figure 2 Solid line depicts the steady-state tenq>eratures corres- 
ponding to the climate solutions for the zero-dimensional 
climate model with variable ice cap as a function of 
solar constant in units of its present value, Tlie Roots I, 

II and III correspond to those in Figure 1. The dashed 
curv,'c shows the solution that would be obtained in the 
case of infinite horizontal heat transport in which the 
planet is isothermal. 

Figure 3 Sine of steady state ice-cap edge latitude x versus solar 
constant in units of its present value for a one-dimensional 
climate model with no horizontal heat transports* cf* 

Eq. 17. 

Figure 4 Temperature Cc) versus sine of latitude for the cases no^ 
transport, infinite transport and earth (schematic). 
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I 

f Figure 5 Mean annual northern hesaspheza aonally-aveisaged teieperature 

I 

I (*C> versus sine of latitude for the two-taodle approximation 

I to the diffusive model (solid line), the observations 

\ (circled dots) end a fit Including F^(x) (dsshed lines) 

* from north (1975b). Cf. Eq. (32) ff. 

Figure 6 Sine of steady-state ice edge latitude a: 'lir 

constant iv mits of its present VAlue for rfhe Budyko 
Jtodel; divergence of iioriznr/cal heat VMz is given 

by (33)© The solution id given b; (34)© 

Figure 7 The Greenes function G^(xvX^), defined by (40)© (41) and 
(44)© for the constant diffusion model versus sine of 
latitude % and for which a heat source is loossted at sine 
of latitude Note that smaller values of diffusion 
coefficient D lead to a more localized responsao 
Figure 8 Same as Figures 3© 5 and 6 except that the model employs 
diffusive heat transport* The model is defined by (15)© 

(22) and <23) with the analytical solution given by (37)* 
Figure 9 Observed surface temperature of the symmetrlZ^:J northern 
hemisphere (dots) and the representation of the surface 
temperatures obtained with the 00, 11 and 20 modes listed 
in Table 1 (solid curv^e). The surface temperatures are 
for northern hemisphere winter and spring, but because the 
temperature fields have been sytaenstrized, the temperatures 
for summer and fall are obtfiined by reversing the absci?=isas 
(North and Coakley, 1979)* 



Figur« 10 (^served i 2 tfrs»«d flusas estitted by ayassttrlsed tu>rtheT*n 

hesisphate (docs) and the representecion obtained with :| ' 

the 00c XI ^d 20 asctdes listed in Teble 1 (solid curve). | 

The dished curve ehoua the fit Obtained by using the 00. ^ ^ / 

11 and 20 leodea of the temperature field (Table 1> in (53) 

(KortU and Coakley, 1979). I 

Figure II Observed albedo of the synzustrlred northern heisiaphere (dote) ' 

end the representaticn of the albedo obtained with the 00. I 

T 

11 and 20 taodea listed la Table 1 (solid curve) (Horth and :j ^ 

Coakley, 1979). 

Figure 12 Distribution of incident solar radiation (dots) and the 

representation of the distribution obtained with the 00, | ; 

11 and 20 ecdes listed in Table 1 (North end Coakley, 1979). 

Figure 13 Teopleths of decrease in global average teaperacure calcu> 

lated for IT decrease in eolar constant (Coakley, 1979). t 

Figure 14 Equillbriusi position of ice line (-lO'C iaothera) as a ! • 

function of solar constant. Results obtained with radiation 

■ f 

panmeterlxatlons adopted by North (1975b) (solid curve) and 

i 

by Coakley (1979) (daahed curve)* For co's^ariaon -10*C 

Isothcra of lowest GCH level from calculations perfotssad /' 

by Wetheraid and Manabe (1975) are represented by points 
(Coakley, 1979). 

Figure 15 Scabillly of northern hemisphere Ice sheet as a function of 
its latitudinal extent for various cos&inations of .bllquity, 
eccentricity and precession (longitude of petihelion 

ssaaaured froia cho position of the northern heolaphore < 
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Hlnt«t Bol8tic«). To tli« right ot the •*s®ro-r«8l»e’* «**« | 

m&u «c.^«m»Ufeion « 3 tc 9 ®ds ablation end th« ice sheet grows i | 

to the S«it eblatlon doalnate?. end the Ice aheot ahrlnVa. | 

8 ^ot»8 • Bt&ble e<iwillbrl«mj 0 an unstable equllibriua ^ j 

0. Pollard, Kature . 271 , 233. 1978. Cop 2 ?righted by ■ :-;| 

* ■ 

»^ac«illa« doumala Ltd., leptinted by P«r®i»olon). . | 

Figure JOe ClobKl average t*»p®ratoro as a funetlon of solar constant | 

■? 

for 5 values of ronslly esymetrlc forcing («»*). | 

Figure 16b Sine of the latitude of the -10"C laothera «« a f^mction 
of longitude for the 5 value* of Myr^netrlc forcing 

(D. Bextaao and D. A. Short, Jj[. £SijL* 519, { 

1979. Copyrighted by Ataerican M^iteorologlcal Society. 

Reprinted by permieaion). 1 ^ 

Figure 17 Ihe potential function for the finite asspllt^de atablUty 
analyelo of the aero-dl®»n8 tonal cllmte i»<feel defined by 
Eqa. (8), <13), (75) and (82). The extrena are labeled 

the usM itt Fi|ture«s I ftnd 2* 

.* • 

Figure 18 Schematic graph r,f Che right-hand aide of E«s. (105) 

(denoted f(l) in the figure) veraua the staMlity parameter 
X. Intereectlons with the (flat) <lQ/dx^ grasph Indicate 

discrete elgenvaluaa of the ayate©. Ntate that If 

dQ/dx^ 18 negative, the lowest eigenvalue Burnt be negative. | 

is^tlylng instability. 

Figure 19 Contours of the potential function (lU) fo«‘ « two-Bode 

one-dltaensional climate «»del are plotted Ire the I^,l 2 ; 

plane, the tvo-»ode em^'lltudes. Ihe atatoa 1. 11 and 111 

correspond to the previous figure*. Fro® Worth ei, al_. (1979) . 


nmr^ 20 Latic«d®-4eF«a4«nt mmitivity (i.e., atatlonery r««pona* 
Jfunctiea) G(a,y) co^uted vith diffusive traoeport, 

8udyko»B Infrared rule and «n iBot'uersssl ice-cap edge 

at m ^0 caaee eho<cri are for perturbations 

at 25,5* latitude (y • 0.4) and at 45* latitude (y • 0.7), 

In both case^ additional heat is absorbed in the Ice edge 
region, and the effect Increases with the global teaperature 
aensltivity (here 8^ • i,6«c p«r IX change in eolar con- 
stant) as well as with the prexltalty of the perturbation 
to the ice edge. 

Figure 21 Schematic illustration of the tluctuatlon dlasipatlon theorem 
in the weak fora. 

(a) a typical time series of global Ceseprature with a ahlfc 
in the mean due to a step change In forcing at t * 0. 

(b) an average of sany aaeple functions siallar to that 

in (a). Averaging allalnatea the fluctuations, Isavlng only 
the mean, which ehsnges froa one stationary value to another 

g tloe fjjg climatic response tiaa, proportional 

to the sensitivity 8 , 
o 

(c) Autocovariance of dT^ - . <T^>. Xhls ftmctlon 

decre,««« over a characteristic lag time the auto- 
correlation tlae for fluctuations in T^. if the fluctuation 
dissipaclcn thcorera holds, t « t , 
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